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CHAPTER I 



INTRODUCTION 



We have attempted to demonstrate the Wigner cusp phenomenon 
in neutron elastic scattering and total cross sections near the thresh- 
olds for excitation of the first excited states of the dominant iso- 
topes of Li, Fe, Zr, Ba, and Ce o The energy region investigated 
extended approximately 300 kev on either side of threshold. 

2 34 5 

Theoretical work (see, for example, Nevton, Baz, ^ Fonda, 

6 , 7 . 

and Meyerhof ; provides detailed predictions of possible cusp effects. 
The effects originate essentially in two fimdamental physical principles 
the conservation of incident flux, expressed in the theory as lonitarity 
of the scattering matrix; and the uncertainty principle of Heisenberg, 
expressed in the theory as analyticity of the scattering matrix. The 
magnitude and shape of cusp anomalijes can yield information useful in 
scattering analysis and in establishing spins and parities of excited 
states . 



Cusp anomalies have been seen in proton elastic scattering near 
thresholds for (p,n) reactions The theory shows that one essen- 
tial requirement for a sharp cusp is that coulomb effects be absent 
from that reaction channel which provides the threshold. Any reaction 
with outgoing neutrons, such as (p,n) (d,n) (a,n) (n,n*), satisfies 
this requirement- 
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The advantage of studying Wigner cusps in elastic neutron^ 
rather than elastic proton scattering, for example, is the easier 
analysis caused by the absence of coulomb effects in the incoming 
channel. After estimating possible cusp magnitudes, ve concluded that 
the effects of inelastic neutron scattering reactions on elastic neu- 
tron scattering might be quite noticeable for selected cases. We 
hoped thereby to find a method for spin assignment of excited states 
and to elucidate interaction mechanisms between nuclei and neutrons 
of few Mev energy. 



CHAPTER II 



THEORT'OF WIGNER CUSPS 

Previous articles on the theory of Wigner cusp phenomena have 
generally emphasized the simple spin cases (e.go, spinless particles); 
also, the results are not alvays given in a form readily compared vith 
experiment. We vould like to supplement the previous results by con- 
sidering more general spin and angular momentum configurations. Our 
results 'will provide a convenient form for experimental comparisons. 

In addition, the form of the results for pure Wigner cusp effects will 
be particularly suitable for subsequent deduction of “energy-averaged*’ 
cusp expressions, as first derived in Ref. 6. 

Pure Wigner Cusp Theory 

We will compute the effect of reactions of the type X(n,n*)X* 

on neutron differential elastic and total cross sections near threshold 

12 

energies. Blatt and Biedenharn have derived necessary general ex- 
pressions for cross sections. For completeness, we give in Appendix A 
a summary of the physical basis for the cross section formulas, start- 
ing with the wave functions and the scattering matrix, 

In applying these fomulas we assume that elastic scattering 
with change of orbital angular momentum i may be neglected. Because 
of the parity rule requiring Ai = 0, 2, 4 . . «, elastic scattering 
with change of & turns out to be an unimportant phenomenon in most 
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k 
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situations. In particular, in the present experiment this restric- 
tion is very mild, since, with one exception, the target nuclei were 
even-even. Hence, with 1=0, 1=^, s=^ (l, i, and s are target, 
projectile, and channel spin, respectively), there is no possibility 
of elastic scattering with change of i. The exceptional case is that 
of Li (l = 3 / 2 ) « In this case, ^ 2 is limited by penetrability 
effects. For examp.le, although incident i = 0 vaves theoretically 
can be scattered elastically into £= 2 outgoing waves (with com- 
pound spin J = l), it happens that in the energy region of interest 
kR = 0.3 any process involving i = 2 is highly improbable. Con- 

sequently, for the present experiment, specializing the derivation to 
= 0 is not a restriction on the validity of the results. 

If one considers angular momenta situations in which there can 
be scattering involving change of channel spin s, but no change in 
channel energy, complexities arise in a derivation of cusp effects. 

To preserve clarity in the presentation, we will restrict our deriva- 
tion to situations in which channel spin s is separately conserved, 
in addition to the previous assumption that orbital angular momentum 
i is conserved. It will be shown that this restriction does not af- 
fect the validity of the derivation as it applies to present experi- 
mental results. We will indicate the nature of the additional complex- 
ities in cusp theory, which are present without this restriction, at 
a suitable point in the derivation. 



The following equations are numbered 3-l6, 4,5^ and 4.6 in 
Ref. 12 o We have specialized them to elastic scattering with no change 
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of a, Z, s or i . Accordingly, the differential elastic scattering 
cross section is given by 



d0 



a; a 



I+i 
= E 
s=|l- 



(2s + 1)/[(2I + l)(2i + 1)]^ da 



as; as 



( 1 ) 



vhere 



da 



OB;as 




( 2 ) 



Bj^(as;as) = (l/U) 




X Z^( Z^ J^, S L) 






( 3 ) 



where a refers to the type of incident particle and the state of the 
struck nucleus; is (2jt) ^ x the deBroglie reduced wavelength of 

the incident neutron/target nucleon system; s, I, and i are channel, 
target, and projectile spins, respectively; and Z, J are orbital and 
total angular momentum, respectively. Subscripts 1 and 2 are dummy 
indices for use in summation. All sms are unrestricted, except that, 
in practice, only one such sum, for example, runs to infinity, be- 

cause of the vanishing of Z coefficients. Re stands for real part 

of the bracketed expression . S"^ . „ is the S-matrix element for 

asZ^asZ 

elastic scattering with Z and s conserved ( see Appendix A) . 
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L-i +i^ i i i i 

Z(ii ^2 ^2^ s L) = i ^ ‘^{2£^ + 1)2 (2ig + l)^ (2J^ + l)2(2J2 + l)^ 

X W(i^ ^2 •^2' ® ^2 ° ^1-^1 -^2 ^ 

l4 15 

where W is the Racah coefficient. Tables of Z are available » 

l6 

The properties of W are discussed by Biedenharn, Blatt^ and Rose. 

The energy dependence of the elastic cross section is contained 
in We will use the unitary and analytic properties of the S-matrix 

to obtain an expression for the elastic cross section in the vicinity 
of a reaction threshold. 

Among the inelastic reactions which can proceed from a given 

incident partial wave, we consider only the one with lowest threshold 

energy. Note, however, that an excitation of a higher excited state 

is not excluded as long as its threshold has the lowest energy of those 

reactions which may proceed from a particular incident partial wave. 

3 

Baz showed that the procedure for including higher inelastic thresh- 
olds of the same dominant partial wave is straight-forward in principle, 
although much more complicated. Since the present experiments were done 
in energy regions near first excited states, the case we consider is 
adequate, and the physical content of the results is not obscured by 
these complications. 

Also, other true inelastic processes (e.g., neutron capture) 
occurring at lower energy than the threshold energy of interest are 
neglected. 

In addition, we restrict the energy region in which the deriva- 
tion is valid to that in which the reaction is dominated by a single 














f, 
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incident partial vave. In practice, such an enei-gy region is usually 
sufficiently large to be of experimental interest. For example, an 
incident d wave contributes more than 90 percent of the cross section 
over an energy interval of 15 O kev above threshold for inelastic neu- 

+ i^to 17 

tron scattering leading to the 2 first excited state of Ce . One 
can explain this in terms of the optical model by noting that Ce has 
a large transmission coefficient for (low energy) outgoing s wave 
neutrons^ and that sufficiently close to threshold^ inelastic neutron 
scattering is alvays dominated by that incident partial vave vhich is 
associated (through angular momentum conservation) vith an outgoing 
s vave. 

Finally^ to further simplify the presentation^ but not as an 
essential limitation on the theory^ ve will only consider reactions 
dominated by outgoing s vaves. The complete theory in this respect 
shovs that reactions dominated near threshold by outgoing vaves of or- 
bital angular momentum greater than zero^ produce cusp effects 

in the (i^)th derivative of the elastic cross section. Therefore^ 
these cases have less experimental interest. With one exception^ this 
restriction is not essential since ve measured threshold effects of 
reactions vith outgoing s vaveSu The exceptional case is the reaction 
Zr^^(n^n*) Zr^^ (l-75 Mev) vhich Tucker^ et have found dominated 

19 

by outgoing p vaves already at about 30 kev above threshold • The 

90 

necessary modifications to the theory in the case of Zr will be ob- 



vious . 
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Under the above assumptions the reaction [see Eq. ( 5 ) below] 
may be described by its partial reaction cross section 



J' 

""ra'i's* 




,J' I 

'’a'i's'ja'i's' 



(4a) 



where primes identify the particular channel quantities vhich couple 
to l' , the dominant incident partial wave and g 1 - (2J' + l)/(2s' + l) 
is a statistical factor obtained after coupling angular momenta, averag- 
ing over initial magnetic quantimi niombers and summing over final magnetic 

I'D 1 p 

quantum numbers (see Blatt and Weisskopf or Blatt and Biedenharn )o 
To simplify the notation, we now drop a number of subscripts. 

First a' is superfluous since we always consider a neutron incident 
upon the groiind state of a target nucleus. Second, we reduce the double 
subscript set on the S-matrix element to a single set; the meaning of 
S, as given in Appendix A, should be kept in mind. 

Equation (^a) becomes: 




Equation (4b) is essentially the statement that the scattering matrix 
^ is uniuary (see Appendix A). The corresponding expression for the 
total reaction cross section is: 



a s a ., = 
r 



I+i 

E 






2s * + 1 ' o jt 

tii.il Isr-rmarvTj 



Inverting Eq. (4b), the S-matrix element may be written; 






(6) 
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J * 

Equation (6) defines 26, i i as the (real) phase of the S- 

Jj s 

J * 

matrix element S.i i . To simplify the notation, we have expressed 
Z s 

J ’ 2 J ' 

a , I in \anits of jr ife r t • This normalization will be used in 
the remaining discussion except in instances where clarity requires 
that the units be specifically restored. 

Equation (6) is certainly valid above threshold. The usual 
procedure is to expand the complete phase (which may have real and 

j * 

imaginary parts) of S , i in energy about threshold. One consequence 

Jo S 

of this procedure is that the fixed value of the real phase at thresh- 
old appears in the final expressions. On. the other hand^ Eq. (6) is 
quite general and there is no reason to assume that the form of Eq. (6) 
should not persist even away from threshold (above threshold). Conse- 
quently^ we shall allow the real part of the phase of i.e.^ 26, 

as well as to vary over the energy region of interest, if the 

physical situation requires that they do so. This leaves us free later 
to examine the consequences of energy averaging over a cusp region where 
the phase 26 and the cross section might be fluctuating rapidly 

because of the presence of many narrow resonances. 

7 

Using R-matrix theory, Meyerhof has shown that the S~matrix 
can always be written in a form equivalent to Eq. (6) if ~ inde- 
pendently whether the energy is above or below threshold. An analytic 
continuation in energy across threshold always exists, even in the 
presence of many narrow and possibly overlapping resonances which would 
produce rapid variation with energy in the phase and in the reaction 



cross section. 
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The prescription for continuing Eq. (6) below threshold, assum- 
ing outgoing £ wave neutrons in the reaction, is given by Eq. (29a) in 
Ref. 7; essentially, one may let = |a^j above threshold go to 
i|a^j below threshold. It is interesting that in R--matrix theory this 
prescription actually results from a change in the effective logarithmic 
derivative of the reaction channel wave function from a pure imaginary 
number above threshold to a pure real number below threshold. This is 
related to the fact that outgoing wave function (exp ik^r)/r has to 
change to (exp ~|k^jr)/r as the energy (E) is changed from above 
threshold (E^) to below. =n/ 2M(E •• E^)/fi^ , M = reduced mass 

of the outgoing neutron. This prescription becomes more complicated if 
one considers outgoing waves in the reaction channel with 0 . The 

necessary modifications for extension of the present theory to these 
cases is also given in Ref. 7* 

Using the fact that i < < 1 near threshold, we can 

1C Xf S ’ 

simplify Eq. (6) to 

and, upon substituting the analyticity prescription into Eq. (7), we 
obtain 



. (exp 2i - i I J ]■) (8) 

2 r A 1 

where^ following Newton^ \ B J ^ ^ ® 

E < E ^ where E is the reaction threshold energy. It should be 



n 
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J ^ 

noted that a , in Eq. (8) is strictly equal to 
r i/ s 



the accompanying 
sion. 



,( |e - E 1 ) 
ri's ' a O^' 

contains the analytic properties of the expres- 



,j' 



We have separated S., i into two terms. The first, 

i/ S 

J ’ I, 

exp 2i 6 . r I ^ can be interpreted as the "S-matrix element -without 
z s 

n 7 

reaction effects," if £ =0. When this term is substituted into 

Eqs. (l), (2), and (3)> and the indicated svims are performed, the re- 
sult is the "differential elastic cross section without reaction ef- 

The second term. 



fects, " which we designate by 



J * J ^ r 1 ^ 

-(exp 2i6.i , . r produces a cusp term after substitution 

i/S ri/S L^J 

into the same equations. 

Accordingly, the Wigner cusp term A is defined by the 

set of equations 



da = da + A da 

a; a a; a a; a 



(9) 



where 



I+i 

s'= I-i 



( 10 ) 



^ ' .rvs ' = /(2s' + 1)] Z AB^(ob';Os') P^(cos 0)dfi ( ll) 



as ;as 



L=0 



4 • 






- m 











Shi 



A 
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AB^(os':as’) = (lA) ESS Z^(i, J, £* j', s’ L) 

'>• Jl il J' ^ ^ 



X Re 



(1 - exp 2i (exp 2i | 



+ S S S Z^(i’ J' i J , s' L) 
J- ^2 



X Re 



(exp 21 6^,^,)* i (1 - exp 21 



J * 2 J * 

vhere a , is still expressed in units of rt g , [see Eq. (6)]< 
Z* X/ s s 

The sums over £^f J^, £^y in Eq. ( 3 ) are reduced because of the 
assumption that a single wave f ’ is dominant in the reaction. 

Because Z(i^ *^ 2 ’ s L) = (-1)^ '^2 '^1 ^ 

Re(A* B) = Re(A B^), we may combine these s\ams^ changing dximmy indices 
and dropping the extraneous index a for simplicity^ to obtain 

A BAs'js’) = (lA) S 2 2 2 Z^(i J £' j', s' L) 

^ J i J' 



X Re 



(1 - exp 2i (exp 2i | j J 



(lA) 2 S S Z^^(i J £’ J', s’ L) 
J i J’ 



X 2 sin d'^^, Re 



i-!^|- exp i (26^, g, 






Tahing the Re [ ] in the last equation yields 
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A Bt(s',s') = - (l/U) ZEE Z (i J J’, s' L) 
^ J i J' 



X a 



|sin(26^,^, 


-1. 


4os(26'J,g , 





( 12 ) 



Within the limitations imposed diiring the derivation, Eqs. (9)^ 
(lO), (ll), and (12) describe expected Wigner cusp effects. The cusp 
effect is completely contained in the coefficients in Eq. (2); 

that is, B = B° + A B , where A B is given by Eq. (12). 

Li Li 1 j Li 

Since it is frequent practice (cf. Ref. 20 ) to analyze neutron 
scattering data in terms of the coefficients, the utility of Eq. 

(12) for cusp interpretation is clear. Should a cusp effect he found 
in only one or two of several experimental coefficients then the 

properties of the Z coefficients severely restrict the participating 
partial waves . 

We now recall that the cusp predictions, Eqs. (9) to (l2) are 
not generally valid if scattering is possible which conserves energy 
but not channel spin, because most detection systems are only energy- 
selective. In such cases, the threshold reaction is coupled (via 
unitarity of the S-matrix) to other off-diagonal (different initial and 
final state channel spin) elements, as well as a diagonal element 
("pure” elastic scattering) of the S-matrix (see Eq» (4a)). An addi- 
tional sum over final state channel spins is required in Eqs. (l) to 
(3)* As a result additional cusp terms appear in Eq. (ll). Unless 
polarization is also measured, comparisons of theory and experiment 



Ik 



are ambiguous in such cases . None of the present experimental results 

required this extension to cusp theory for comparison purposes. 

The expected Wigner cusp in total cross sections (a^) is 

easily obtained if ve recall that, by assumption, the interaction in 
th 

the (^') incident channel is purely elastic below threshold. 



f (da/dn)dfi + 






/ (A da/dfl)dft + 






(13) 



where is the total cross section "without reaction effects." 

Using the fact that / P (cos 6)di2 = 6^ (where 6^ 

h J_j, 0 L, 0 

is a Kronecker delta) we obtain 



(2s* + 1) . fl 

' s' = |l-i| (21 + l)(2i + 1) (2s* + 1) ^ V® ^ ^ ‘"ri'lq 



(14) 



Substituting Z(i J £' J', s' O) = 6^^, &jj,(-l)'^“® (2J + l)^ , where 

K K 

6^^, and &jji Kronecker deltas, into the equation for A 

(with L = O), Eq. (12), and putting this resTilt into Eq. (l4), one 
obtains 



Aa - (2s' +1) A»c^(2J'+1) J' 

' s'=|l-i| (21 + l)(2i + 1) j.f|^._s'|V (2s'+l) ‘^ri's' 



- 2 sin^ ^ 

J, J °ri'l0 

- sin 26^,^, 



( 15 ) 



15 



J * 

where we have explicitly indicated the units of (see Eqc (6))o 



Finally^ recognizing that the summations in EqSo ( 15 ) and ( 5 ) are 
identical, we may combine 



Equation ( 16 ), combined with Eq, (13)^ is the general expression for 
Wigner cusp effect in the total cross section, within the limitations 
imposed by the derivation. 

Wigner Cusps — Spinless Particles 

We now extract the special case of spinless particles from these 
expressions, for comparison with other work« 

If all spins are zero, including that of the residual nuclei 
in the reaction, then Eq» (12) becomes 




and obtain the result 




V ^^( 2 J' + 1 ) J' 
(2s ' + 1) ' & 







( 16 ) 



A = "(1/4) Z Z^(i 0 0, 0 L) 2 sin 6 ^ 



t 



s±n(.2&o - 



os(28q - 5^ 
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Since Z(£ i 0 0, 0 L) = (2i + l)^ (i 0 0 0 ] i 0 L O) 



K K 

= , where 5._ is a Kronecker delta, 

£Li X/\a 



the sum in A is reduced to one term, and using Eqs. (lO) and (ll) 

Li 

ve obtain 



00 r 

A da = - (>c^A) 2 a 2 sin 6 \ 

£=0 ^ ^ ^ 



sin(26Q - 6^) 



cos(25q - 6^) 



} 






where 6^ = 5 ^q- This expression can he written as 



A da = 2 Re 



(a° n 

(2 1 ^ (1 - exp 2i S^)P^) 2fl % 2 ^0 1i 



thus 



A da = 2 Re 



UA) f*(e) I exp 21 Sg |J-| 



(17) 



where f*(6) is the complex conjijgate of the scattered amplitude, with 

its explicit form evident by comparison with the preceding line^ and 

2 0 

where (a ) = (it k i.e., the experimental reaction cross 

^ r exp rO 

section. 

U 

This is just the cusp term obtained by Baz, for the case of 
spinless particles. Baz emphasized the explicit appearance of the am- 
plitude, suggesting that there is additional theoretical significance 
in cusp data beyond cross section measurements (squared amplitudes) 
that are usually made. We woxild like to caution that most cases in 
practice involve non-zero spins; consequently scattered amplitudes, 
in general, are azimuth angle (9) dependent. Comparison of theoretical 
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cusp effects with experiment may be undertaken if the theoretical ex= 
pressions have been averaged and summed over magnetic quantum numbers. 
This is essentially what we have accomplished, and, in general, little 
vesitge of complete amplitude forms is found in the resulting expres- 
sions . 

On the other hand, it is useful to see from Eq. (17) that one 
can estimate expected Wigner cusp magnitudes from the e3<pression 



the remainder of the cusp expression being of order unity. The ex- 

1/2 

pected percentage cusp effect is inversely proportional to (da) ' 
and exploratory experiments might well he mdertaken at angles where 
the cross section is leasts 

The spinless particle result for total cross sections is a 
trivial reduction of Eq<> (l6) to 



1 



|a do| = o^(da)2/x 





This concludes our discussion of p\ire Wigner cusp phenomena. 



Two other special cases, applicable to the present experiment, are de- 
duced in Appendix B. 
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Theory of ^’Energy- Averaged*’ Wigner Cusps 
We turn now to the problem of calculating the expected threshold 



effect, if the measured cross sections are effectively energy averages 
of the actual detailed cross sections, containing closely spaced reson- 
ance structure. That the resulting effects are related intimately to 



showed that the energy average of Wigner cusp expressions leads to a 
predicted threshold effect "vdiich we shall call an ’’energy-averaged" 
Wigner cusp. 



and scattering cross sections are dominated by many closely spaced 
resonances. By closely spaced we mean a large mamber within the energy 
interval A E over which the averaging is performed: either experi- 

mentally, because of finite incident beam resolution or because of 

overlap of the resonances; or theoretically, by choosing A E > > D, 

7 

where D is the average level spacing. No essential restrictions 
will be placed on the level width to spacing ratio, p/d. 



effects is contained in ABj^(s*;s*) (see Eq. (12)) we can isolate the 
energy dependent quantities: 




The physical situation we are considering is that the reaction 



Noting that the energy dependence of differential elastic cusp 




J 



J 



( 




Brackets, ( ), imply that an energy average is to be taken 
over an interval containing many resonances. The phases and 

X/S 



19 



J » 

5 - s I are rapidly varying with energy over the averaging interval and 
Z s 

we will assume they are uncorrelated for J ^ J* and/or Z f This 

is the well known statistical assumption for the compoiind nucleus* In- 
tuition might suggest that uncorrelated phases would lead immediately 
to breakup of expression (l8) into a simple form, for example, the 
product of the averages of the factors* It turns out that such a break- 

j * 

up is by no means obvious, because of the presence of cr ,, , within 

r i/ s 

the averaging brackets. 

Using R-matrix theory and applying the random sign assumption 

7 

to certain value quantities y > which will be discussed, Meyerhof 

K, 

showed that 





. 2 

sin 

in 26 





J’ 

^ri's' 

0 



1 



(19) 



and 




- ^ fo\ 

• orJ' j ^ ~ loj 

Sin 26^. s. 



(20) 



The value quantities y mentioned above were introduced by 
21 

Wigner and Eisenbud to parametrize the terms of the R-matrix, each 
of which has a resonance denominator. Thus, these quantities corre- 
spond to resonances and could be indexed in our notation by > 

where X specifies a particular resonance of the set associated with 
angular momenta J, £ and spin s. The random sign assumption leading 
to Eqs. ( 19 ) and (20) is made for value quantities of different index 
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J, £, or s, but. is not justified for relating value quantities of 

the same J, i and s, hut different X; i.e., different resonances 

initiated hy the same partial wave may have correlated phases. In 

particular, the reaction cross section appearing in the averaging 

brackets in Eqs. (19) and (20) proceeds through resonances correlated 

in phase to the remaining expression in the bracket, and it is surpris- 

Y 

ing that the simple results, Eqs. (19) and (20) were obtained. 

In Appendix C, we show how Eqs. (19) and (20) can be applied 
to expression (l8) to give the result 



where 6 is a Kronecker delta. 

Part of the proof of Eqs. (19) and (20), which led to Eq. (21), 
depends upon a cancellation of path integrals used in the energy averag- 
ing process, which will be complete only if normally slowly varying 
quantities remain slowly varying throughout the averaging interval. 

For the case of outgoing ^ waves in the reaction channel, one of these 
quantities (the logarithmic derivative of R-matrix theory) does not 
vary slowly right at reaction threshold. In essence, this difficulty 
is just the result of rediscovering pure Wigner cusp phenomena as one 
formally approaches very close to threshold in the averaging process. 





( 21 ) 
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This means that the derivation of Eqs. (19) and (20) is not strictly 
valid arbitrarily close to threshold. The actual experimental result 
would depend on the exact energy distribution of the beam. As long as 
the energy distribution of the beam does not overlap the threshold, 
the results, i.e., Eqs. (19) and (20) and consequently Eq. (2l), will 
be applicable. 

After justifying Eq. (2l) we have completed the derivation of 
''energy-averaged" cusp expressions. Recalling Eq. (12), which contains 
expression (l8), and using Eq. (2l), we obtain 

(^J' 

Z J’ i' J’, s' L) I 

J’ 0 



( AB (9';s') ) = - (lA) 




the sums over i and J have vanished because of and 5jj» 

in Eq. (2l). The "energy- averaged" differential cross section cusp term 
is obtained by substituting Eq. (22) in Eqs. (9-H). 

Substitution of relation (20) in Eq. (l6) gives the expected 
result for the total cross section. 







(23) 



That is, there is no "energy- averaged" cusp in the total cross section. 
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This is consistent with the basic optical model assumption, used in 
relating the model to experiment, that the energy- averaged total cross 
section, being linearly dependent upon the scattering matrix elements, 
is smooth in energy even across reaction thresholds. We will come back 
to this point when summarizing the theory. 
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"Energy-Averaged" Cusps — Spinless Particles 

As an example, for comparisons to previous work, consider the 
special case of spinless bombarding and target particles 

l=i=s'=0; i'=J' 



Equation (22) becomes 



2 

( A Bj^(0;0) ) = - (1/4) Z {£' V V i', 0 L) -j 



^ri'O 

0 



(22^' 



Using Z^(i' i' Z' Z' , 0 L) = (2i' + l)^ (Z' i' 0 0 ] Z' Z' L O)^, 

indicating explicitly the mits of \b «i^)^ S’lid l>y substituting Eq. 

r i/ u 

(22s) into Eqs. (9-ll), we obtain 



( Ad 0 ^ ) = - (2Z' + 1)(« k^(2Z' + 1) (a-^' )/4n)Z (Z’ Z' 0 0 \ Z' Z' L of 

CCj Ou L 

X d2 

^ - (2r . (2M 

where (0 „i) means units have been restored to (0 „,) and experi- 
rZ exp rZ ^ 

mental values may be substituted for this quantity. Equation (24) is 
the same as the formula, Eq. ( 16 ) of Ref. 6 . 

"Energy-Averaged" Cusps — Spin Incident Particles 

We give the results for a second special case of Eq. (22) of 
particular interest in comparison with the present experimental data, 
where the incident particle is a neutron, the target nucleus is even- 
even, and incident d waves dominate the reaction. 
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I = 0, i = -I, J’ = i' + -I, i' = 2, and s’ = ^ 



i'+ i 



< Ado(r -2)> -- ,, i®j' 

J —‘Jj "2 



X £ [Z^(2J' 2J', I L)/(2J' + l)]Pj^|Q|dft (25) 

T I 

where R^, = (a i) /(cr ,) 

J r ^2 exp' tZ exp 

Although not immediately obvious, Eq. (25) is identical to Eq. 

(l6') in Ref. 6. To compare, one must identify R»i. i. = ( i '+!)/( '+1 

i/ ■*" 2 

+ i'r), R.i i=i'r/(i'-+l + i'r), where r is defined and used 
i! - 2 

in Ref. 6. 

Since Eq. (25) is important to present experimental comparisons, 

it will be useful to give the terms in detail. It turns out that the 

cusp prediction, Eq. ( 25 ), is quite insensitive to the ratio R , ; 

recall that R , was defined as the fraction of the reaction probabil- 
J 

ity vhich originates from compound resonances of a particular J', for 

example t' + We make the reasonable choice R , = ^(2J’+l)/(2i'+l) 

J 

vhich is the same as the fraction of the total number of incident 

th 

states of the ( i ' ) partial wave which have total angular momentum 

J'. With this choice of R^i, Eq. ( 25 ) becomes 

J 



(a d(j(i’ =2)) = - ((0 ^ 

TX, 




•> 1 

" 2 



“ ■ {^0} 

X [1 + 38/35 P 2 (cos 0) + 18/35 Pj^(cos 0)]dfl 



(26) 



2k 



Equation (26) actually applies to any scattering experiments 
on medium and high atomic numbered even-even targets at energies where 
random phase concepts are useful^ and near the threshold of a reaction 
dominated by incident d waves « The outgoing particle in the reaction 
channel must be a neutron, as mentioned previously, in order to obtain 
a sharp cusp. 

In conclusion, we remark that the expressions for "energy- 

averatged" Wigner cusps depend only upon statistical assumptions, and 

not upon a particular model, for example, the optical model; rather, 

these expressions simply give the change in cross sections expected at 

reaction thresholds, due to unitarity of the scattering matrix. Meyerhof 

6 

showed, however, that these statistical assumptions are consistent with 

the optical model assumption that the scattering matrix element consists 

22 

of a smooth plus fluctuating part. 

Also, we note that "energy-averaged" Wigner cusp expressions 
predict no below- threshold effect. Consequently, in regions of mass 
number and energy where "energy-averaged" theory applies, it is not 
necessary to specify the dominant outgoing wave in the reaction channel, 
in contrast to pure Wigner cusp theory (see the discussion just before 
Eq« (?))• In particular, outgoing 2 'wave dominance of the reaction in 

90 

the Zr experiment, discussed earlier, causes no difficulty in analysis. 
Determination of the dominant incident partial wave in the reaction and 
Eq. (26) suffice. 



CHAPTER III 



EXPERIMENT 

Method 

Neutron Differential Elastic Cross 
Section Measurements 

Neutrons were produced by the T(p,n) reaction^ using tritium 
loaded zirconium and titanium targets for the proton beam from a High 
Voltage Engineering Company 3-Mev Van de Graaff generator stabilized 
to about 1 kev. Various tritiimi targets^ described in Table I^ of from 
10 to 70 kev thi.ckness at 1.02 Mev proton energy (as determined by the 
rise method at threshold) vere used. An eccentric circular motion was 
imposed upon the tritium target. Consequently^ the proton beam struck 
a ring-shaped area on the target^ which minimized cooling difficulties 
and also effectively averaged over nonimiform tritium deposition in 
the target^ thereby producing a time integrated neutron flux very 
nearly proportional to accumulated proton charge. 

TLie neutron beam was collimated by a rectangular conical open- 
ing in LiC 02 *“loQ'^^^ paraffin source shielding of characteristic thick- 
ness 18 in. The geometry and shieldi.ng of the experiment is shown in 
Fig. 1. The collimation angle was such that the energy spread in the 
neutron beam due to T(p,n) kinematics was approximately k kev. A 
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standard long counter ^ placed at 0® with respect to the beam (see Fig. 
1 ) was used to monitor the neutron yield. 
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Fig. 1. Experimental Arrangement — Plan Views 1. Tritium 
loaded target; 2. Air-water spray coolant; 3« Eccentric wobble for 
additional cooling effect; 4. Paraffin shielding (loaded with 
LiCO^); 5* Paraffin (LiCO^ loaded) rectangular conical collimator; 
6. Scatterer in "in Beam" position; 7* Scatterer in "Out" position; 
8. Anthracene or stilbene crystal (a right circular cylinder, 2" 
dia X l" long); 9» 6810A RCA photomultiplier tube; 10. Preampli- 
fier (linear and nonlinear gamma suppression outputs); 11. Long 
Counter (BF^) in monitor position. 
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The neutron beam was allowed to strike the scatterer^ which was 
made slightly smaller than the aperture of the neutron beam, at about 
10 in. from the collimator end. The lithiijm scatterer was a right cir- 
cular cylinder oriented vertically; the remaining scatterers used were 
in the shape of a slab, oriented with their normals parallel to the 
floor and at 40^ to the incident beam direction o The scatterers used 
are described in detail in Table II. Scatterer thicknesses were be- 
tween 1/5 and 5/8 neutron mean free paths at the energies of experi- 
mental interest. These rather thick scatterers were necessary to obtain 
sufficient coiinting rates . Scattered neutrons were detected at two 
angles by organic scintillation crystals (2 ino dia x 1 in. long anthra- 
cene and 2 in. dia x ^/k in<, long stilbene) placed one on each of the 
transmission and reflection sides of the scatterer at distances varying 
from 6 to 10 in. and with the crystals oriented as shown in Fig. 1. 

The scatterer to crystal distance varied with detection angle and was 
set at the maximum consistent with a neutron background/effect ratio 
less than two« A saturated space-charge pulse-shape discrimination 
output of the detection preamplifiers was used to suppress gamma-ray 
background. Typical neutron pulse spectra, given in Fig. 2, illustrate 
that we also used an integral discriminator to bias out inelastic neu- 
trons . 

To facilitate neutron backgro\ind subtraction, the scatterers 
were suspended on thin wires from an overhead crane so that the scat- 
terer could be alternately removed and replaced in the beam at each 
neutron energy (see Fig. l)<. 
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TABLE I 

Description of Tritixan- Loaded Targets 



Index 


Made by 


Yield 


Thickness 






neutrons/ sterad 
microcoulomb 
2.50 Mev 


(at T(p,n) thresh- 
old 1.019 Mev) 
kev 


1 


HAR 


1.5 X 10^ 


22 


2 


HAR 


0.37 X 10 ^ 


7 


3 


HAR 


l<-.8 X 10 ^ 


60 


k 


OR 


0.33 X 10^ 


11 


5 


OR 


1.1 X 10 ^ 


11 



HAR = Harwell, England (Tritium-loaded titanium) 

OR = Oak Ridge Natl. Lab., Tenn. (Tritium- loaded zirconium) 



TABLE II 

Description of Scatterers 



Scatterers Used in da Measurements 



Material 


Dimensions 

(inches) 


Chemical 

Purity 

ii>) 


Angle between 
slab and beam 
( degrees ) 


Cerium a. 


k ,0 X 3.0 X 0.63 


99.8 


50 


b. 


3.25 X 2.75 X 0.19 


99.8 


50 


Barium 


k.o X 3.0 X 1.19 


97.5 


45 


Zirconium 


4.0 X 1.5 X 0.56 


98 


50 


Iron 


4.0 X 2.5 X 0.52 


98 


50 


Lithium 


3.05 dia X 4.4 high 


~ 100 


Rt. cir. cylinder 


Cerium 

Barium 

Iron 


Samples Used in a^( Transmission) Measurements 
1.0 X 1.25 X 0.935 long 

1.0 dia. X 1.87 long dimen- 

Sion was parallel 

1.0 X 1.13 X 0.97 long to the beam axis. 
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Fig . 2 . Pulse spectra from detection equipment at neutron 
energy of O .5 Mev. a. Gamma-Ray-Suppression Spectra (nonlinear pulses): 
1. With the detector (2*‘ dia X 3A” long stilbene crystal) in the neu- 
tron beam and with the pulse height analyzer (PHA) recording all preampli- 
fier nonlinear output pulses. 2. With the. proton beam blocked by a beam 

22 

stopper from hitting the tritium target, with a Na source taped on the 

crystal, and with the PHA recording all preamplifier nonlinear output 

pulses. 3« With the proton beam on the tritium target, with the detector 

22 

crystal in the neutron beam, and with a Na source on the crystal. Two 
integral discriminator biases were set, one on each of the linear and 
nonlinear spectra, and the PHA had been adjusted to accept nonlinear 
pulses which were associated with a time coincidence of discriminated 
linear and nonlinear pulses, b. Linear Spectra (linear preamplifier 
output pulses associated with a time coincidence of discriminated linear 
and nonlinear pulses): A. With the detector crystal in a typical experi- 

mental position ( 90 ® lab) and the scatterer in the experimental position 
in the neutron beam. B. With the detector in the same experimental posi- 
tion as in A, and with the scatterer retracted to its *’Out" position 
(see Fig. l) . C. (A - B) = elastically scattered neutron pulse spectrum. 
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PULSE HEIGHT RELATIVE UNITS 



32 



Neutron Total Cross Section Measurements 

The same neutron beam, source shielding, collimation, and gamma 
suppression described above were used in measuring total cross sections. 

Neutrons were detected in this case by a stilbene crystal with 

dimensions l/2 in. dia X l/2 in. long mounted at 0® to the incident 

beam. Transmission samples were right parallelepipeds whose area normal 

to the beam 1 sq. in-) was sufficient to shield the crystal from 

direct neutrons. (Neutrons were produced over a proton beam spot area 

of about l/4 in. X l/4 in.) The samples were 10 in. from the collima- 

tion (28 in. from the neutron source) and 20 in. from the crystal. 

Sample thicknesses were chosen to allow approximately 60 percent neutron 

24 

transmission at energies of interest in the experiment. The sample 
and a 9 in. long 1 in. dia. tungsten shield were mounted on thin sup- 
porting rods above the ends of two arms of a rotatable horizontal cross. 
A third arm had only a thin supporting rod attached. Thus it was pos- 
sible to coimt transmitted, background, and direct neutrons in sequence 
at each neutron energy by rotating the cross . 

Apparatus 

It remains to describe the electronics systems necessary to re- 
duce crystal scintillation events to (elastically scattered) neutron 
counts. A detailed discussion of our gamma suppression circuit is con- 
tained in Appendix D. The block diagram. Fig. 3^ will be useful in 
following the present discussion. 
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Fig. 3. 



Block diagram of electronics, showing the function of 



each major component. 
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For each of two detection systems an RCA 68IOA photomultiplier 
and preamplifier (Fig. 22) produced linear and nonlinear (gamma sup- 
pression) outputs in response to a crystal scintillation o These two 
outputs fed into a double delay line and RC coupled linear amplifier, 
resp'ectively . 

Each of the amplifier outputs drove an adjustable bias circuit 
(integral discriminator) which generated a standard square 20 volt 
output for each input pulse that exceeded the bias level. After time 
synchronization was established by the use of a suitable delay cable, 
these pulses were allowed to enter a standard coincidence unit. A co- 
incidence output pulse was obtained for each simultaneous occurrence 
of input pulses. Thus, with properly adjusted bias levels on each of 
the linear and nonlinear pulses, a coincidence output pulse signalled 
the detection of an elastically scattered (linear pulse above linear 
bias) neutron (nonlinear pulse above nonlinear bias). These pulses 
were the experimental coxints. 

A 256-channel pulse-height analyzer (PHA) was a valuable tool 
in setting up the gamma suppression circirlt, the amplifier gains, and 
the discriminator biases. We also used the PHA as a monitor on the 
overall electronics performance during a r\m; these monitor spectra had 
a specific role in data analysis for some of the experiments (to be dis- 
cussed). Typical pulse spectra are shown in Fig. 2. The coincidence 
output pulses also served to gate the PHA when spectra such as those 
shown in Fig. 2 were taien. 
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A standard Hansen and McKibben long counter was employed to 
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monitor the neutron yield. Direct neutrons were moderated by a paraffin 

10 3 

cylinder and some were captured by boron in a one-in. dia. BF detec- 
tion tube. The subsequently emitted alpha particle initiated a pulse 

which was amplified and counted as a neutron event. The efficiency of 

-3 

the comter was approximately 2 x 10 and was fairly constant for a 
range of neutron energies from O.5 Mev to 3 Mev. A 1/8 in. cadmium and 
3 in. paraffin shield on the sides and back of the cotmter reduced 
floor scattered neutrons to less than I5 percent of direct neutrons. 
Although absolute neutron yields of our tritium targets were measured 
using the long counter (see Table l), these yields were not required 
in the data analysis. 



Procedxire 

Differential Elastic Cross Section 

Since the experimental interest centered on effects near the 
energies of first excited states, initial adjustments of the apparatus 
were made at these energies. For each experiment the nonlinear preamp- 
lifier controls (see Appendix D) were adjusted for optimum gamma-ray 
suppression, while the crystal was located in the scatterer "in beam" 
position, later called the "calibration position." Care was taken to 
insixre that the crystal plane was parallel to the incident beam axis. 
This aspect corresponded to the aspect of the crystal toward scattered 
neutrons when it was in experimental position (see Fig. l) . After 
these adjustments were completed, the neutron yield and overall effi- 
ciency of the crystal/electronics system were calibrated by counting 
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direct target neutrons at several neutron energies above and below 
threshold. This was done for each of two crystal detector systems. 

After the crystals were moved to the experimental position 
(Fig. l), the same energy range was covered in energy increments of 
approximately one-half incident neutron-energy spread. At each energy, 
counts were recorded both with the scatterer in and with the scatterer 
out of the beam, thus providing a basis for subtraction of neutron 
background. On a succeeding run down in energy the same increments 
were used, but the points interleaved those from the run up; thus, we 
checked data reproducibility. 

All counts were normalized to a fixed amount of accvonulated 
proton charge. When using 30-50 microamperes proton current, we found 
that the necessary counting time was typically 10 minutes for each 
energy point. 

Long coimter monitor neutron, coimts were also recorded at each 
point. Monitor neutron coimts with the scatterer out of the beam 
served both as a check on the condition of the tritiim target, which 
was subjected at times to as much as 70 watts of heat flow, and as the 
basis for a small correction of the data analysis for target yield 
variation . 



Total Cross Sections 

As stated above, the gamma suppression circuit was adjusted at 
a neutron energy corresponding to an excited state of interest. After 
carefully aligning the sample and shield between the proton beam spot 
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and crystal detector, we took a sequence of counts at each energy point 
By rotating the supporting horizontal cross on which sample and shield 
were mounted we were able to place first the sample, then the shield, 
and last the duplicate sample support (for a count of direct neutrons) 
into the beam. After each such set of counts the neutron beam energy 
was increased by approximately one-half incident beam-energy spread 
until a desired energy range was covered below and above threshold. A 
subsequent run down in energy, with interleaving points, verified re- 
producibility. 



CHAPTER IV 



DATA REDUCTION AND ANALYSIS 
Data Normalization 

The basic normalization used in the experiment was to a fixed 
accvimulated proton charge, which implies an assixnption that the result- 
ing time-integrated neutron flux was also a fixed quantity, except for 
the smooth change in the T(p,n) reaction cross section vith energy, 
and except for a ±1 percent uncertainty introduced by our proton cur- 
rent integrator. As mentioned before, the eccentric circular motion 
of the tritium target minimized any imperfection in this assumption, 
by averaging over nonuniformities of tritium deposition. Hovever, 
small variations in integrated neutron fl\ix from the true smooth shape 
of the T(p,n) reaction cross section might have been present vhen the 
neutron energy vas changed, because of changes in proton beam intensity 
distribution vhich occur vhen one sets the accelerator to a new energy. 
If corrections were not made this could have produced scatter in the 
data, as plotted versus energy, in excess of statistical fluctuations. 
For this reason, monitor neutron comts, which had very small statis- 
tical errors, were used to renormalize the differential cross section 
effect counts to a smooth curve drawn through monitor neutron data; 
this curve was assimied to have the true smooth shape of the neutron 
yield versus neutron energy. The normalization was then a percentage 
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correction to each data point depending upon the percent deviation of 
the monitor counts from the smooth curve- This correction was usually 
less than one percent and never more than three percent- It was gen- 
erally impossible to distinguish by eye the difference between corrected 
and uncorrected data, but in seeking cusp effects, we wished to take 
all precautions necessary to establish the relative magnitude of the 
data points plotted versus neutron energy, 

Data Averaging 

The original data for experiments on zirconiiom, cerium and 
barium showed reproducible fluctuations in addition to a well defined 
drop in magnitude over a broad region above reaction thresholds . It 
turned out that energy averaging the original data did not mask the 
essential structure or features. Figure 4 illustrates this point, 
where we have used differential elastic cross section data from a 
ceriimi experiment. Therefore, we were able to improve the statistical 
errors without significant loss of structural detail and also to estab- 
lish a standard energy averaging interval (A E = 50 kev) for presenting 
measurements at many angles . 

Insofar as we were concerned with pure Wigner cusp effects, we 
could not use data averaging procedures in the experiments on lithium 
and iron- Here, the magnitude of the structure in the cross sections 
turned out to exceed that expected of cusp effects; in the experiments 
on iron, the structure was particularly marked- A successful search 
for a cusp effect would depend upon distinguishing this basic structure 



Fig. k. Neutron differential elastic (90“) cross section of 
natural cerium for three widths of energy averaging interval. 

(a) Original data; A E = 7 kev; statistical errors ±5?^. (h) Same data 

averaged over A E = 25 kev; statistical errors 2.5^^. (c) Result of a 

larger averaging interval, A E = 50 kev; statistical errors approxi- 
mately 1^. The averaging in (c) was selected as the standard for cusp 
analysis of cerium. The loss of structural detail upon successive 



averaging is apparent. 
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in the cross sections from a true cusp effect. We illustrate with 
Fig. 5, which gives iron total cross section measurements, and which 
shows that use of 25 kev incident neutron energy resolution washed 
out significant structure elucidated by experiments with 7 kev spread 
in incident neutron energy. 

In addition to washing out structiire, a poor-resolution inci- 
dent beam also has the well known effect of distorting sharp features 
in cross sections. We observed this in differential cross section 
measurements; a comparison of Fig. 11 and Fig. l4 shows a featiire shift 
of about +10 kev as the result of using 30 -kev compared to 7-kev energy 
resolution. This amount of feature shift would make a cusp analysis 
impossible . 

In reporting some of the lithi-um differential cross sections 
we have energy averaged the data for clarity in presentation (see Fig. 
lO) . It is clear from the foregoing discussion that in so doing we 
conceded that there was no cusp effect in the data and that remaining 
structure must be attributed strictly to the elastic scattering cross 
section of the lithium nucleus . 

Data Analysis 

Differential Cross Sections 

We give our calculation of the measured differential elastic 
cross sections using (l) the basic data, which consisted of the counts 
C = Cj '(scatterer in) - (scatterer out), (2) geometrical factors, 
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Fig. 5“ 'Hie total neutron cross section of natural iron, using 
neutron beams of two energy spreads: (a) A E = 25 kevj (b) A E = 7 
kev. The size of the symbol indicates the size of the statistical 



errors . 
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and (3) an overall efficiency measured in a calibration run. In order 
to make a detailed comparison of gross structure, cusp effects in 
cross sections, or angular distributions with theory, or with other 
experimental work, additional corrections for multiple scattering and 
finite scatterer effects had to be applied. These corrections are dis- 
cussed in Appendix E and tabulated in Table III. 

Measured cross sections will be given in the center of mass 
(cm) frame and will be expressed as functions of mean polar CM angle 
( 0 ^(CM)) and mean CM incident neutron energy (E (CM)). Our reason 
for presenting CM data is that theoretical formulas for cusp phenomena 
and threshold energies are more conveniently expressed in CM terms. 

The quantities da(CM), 0 (CM), and E(CM) were computed from the cor- 
responding quantities expressed in the laboratory frame, using well 
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known transformation formulas. In the present experiment the differ- 
ences between lab (laboratory) and CM quantities were negligible except 
in the cases of iron and lithium. 

E(lab) was determined by subtracting one-half the incident 
neutron energy spread from the maximum lab energy of the neutron beam. 
Incident neutron energy spread was computed, assuming statistical in- 
dependence of the two contributions, from (l) the geometrical spread, 
fixed by the neutron collimation angle and T(p,n) kinematics to be 
4 kev, at 1.5-Mev neutron energy; and ( 2 ) the proton energy spread 
caused by the thickness of the metal target, containing tritium. The 
latter spread was computed from the target thickness at T(p,n) thresh- 
old [1.019-Mev proton energy (lab)] as determined from the slope of 
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the rise at threshold and conversion factors relating proton stopping 

power of the target at threshold to the stopping power at the energy 
25 

of interest. Neutron energy spread is directly related to proton 

/ \ 26 

energy spread by TCp,n) kinematics. 

Maximxim incident neutron lab energy present in the incident 
beam was determined from maximum proton lab energy by the same kine- 
matical relations, where the proton energy was calibrated in two ways. 
First, by using the acciirately known T(p,n) threshold energy we ef- 
fectively calibrated the Van de Graaf machine analyzing magnet a,t 
1.019 Mev. The magnet field was measvired with a proton magnetic res- 
onance probe, between the pole faces, to obtain with high precision 
an RF frequency which coiald be associated with 1.019 Mev proton energy. 
We used the intersection of the straight line extrapolation of the 
lead slope of the T(p,n) yield versus frequency c\urve and the frequency 
axis as indicating the trjie threshold. Since proton resonance fre- 
quency and the magnetic field (b) are linearly related (f = p • B/h, 
where f is frequency, p is the proton magnetic moment, and h is 
Planck's constant), and the analyzing magnet is a momentum selector, 

then maximm proton lab energy is E^^(lab) = I.OI 9 f^/f^^ Mev where 
2 

frnrr is the previously described threshold proton resonance frequency. 
In 

12 ^ 

Second, because C has a sharp resonance at 2.07o-Mev lab 

neutron energy, a long counter of the type used in this experiment has 

an anomaly in its efficiency versus energy cxirve at 2,076 Mev due to 

12 23 

the increased cross section of C in the paraffin moderator. Thus, 
a dip occurs in monitor neutron counts at this neutron energy. The 
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proton resonance frequency at this dip established a calibration point 
for the meein neutron energy which served as a check on calibration by 
the first method. The two calibration points conveniently straddled 
the proton energy region involved in the present experiments. 

Initially, the "effect coiants" C were normalized as described 
in the first section of this chapter, giving normalized effect counts, 

^N* 

For a given scattering sample and incident neutron lab energy 
E the relation between Cjj and the differential cross section is; 

C„ = Y f [(E - E„), R ] f JE - E^) n a 
N X R ^ X D R X s 

X [1 - exp(- nff^ d)] da(0, E)/ci^ (2?) 

where Y is the neutron yield of the target per imit solid angle; 
f^[(E - is the intrinsic crystal efficiency at the lab energy 

of the scattered neutron and at the crystal- to-scatterer distance, R^; 
E is the lab recoil energy of the residual nucleus; f_(E - E ) is 
a discriminator correction to be discussed; and are the solid 

angles subtended by the crystal to the scatterer center and by the 
scatterer to the proton beam spot, respectively; nd is the number of 
scattering nuclei per unit area in the direction of the incident beam; 

is the total cross section of the scatterer material; and da(0, E) 
is the differential elastic cross section per unit solid angle. Since 
da(0, E) has been assumed constant [otherwise Eq. (27) would contain 
integrals over solid angle of the finite detector and energy spread of 
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the beam] over the angular region described by and the energy 

region included by the energy spread of the incident neutrons, the 
symbols 0 and E refer to mean polar angle and mean neutron beam 
energy respectively [i.e., 0^(lab) and E(lab)]. 

The corresponding formula for "calibration counts" is 

C = Y f (E,R ) f JE), (28) 

c c X c D 

where is the so3-id angle subtended by the crystal to the proton 

beam spot and is the associated scatterer to tritium target dis- 

tance o 

Taking the ratio Eq, (27)/Eq. (28) yields 

f (E,R ) fjE) fl 

“ f [(E - E_), R j f_(E - E_) n a 
X R X D R X s 

X 0^[1 - exp(- n 0^ d)]'^ 

■where K is a constant normalizing the effect and calibration counts 
to the same accum\ilated proton charge. 

The factors f^^ are necessary corrections vhich arise because 
the discriminator bias on detected linear pulses ■was fixed. As a re- 
sult, the counted fraction of the total detected elastic neutron spec- 
trum is a function of energy. This is sho^wn in detail by Fig. 6 and 
is explained fiirther in the figure caption. In the detection of sec- 
tored neutrons, the factor f^, becomes C(E - E„) because the lab 
— - • JJ D K 

energy of the scattered neutrons rather than that of the incident neu- 
trons determines the discriminator correction. In all except the 
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Fig. 6. Linear pixlse spectra for elastically scattered neutrons 
at three energies, shovring the energy dependence of detection efficiency 
due to fixed integral discrimination bias. The spectra were obtained 
with the detector at 125 “ (lal>) "to a lithiimi scatterer and are the re- 
sult of subtracting the spectrum with scatterer out of the neutron beam 
from the spectrum with scatterer in the beam. This figure only illus- 
trates an experimental point in the data analysis. Actual rvms were 
made with the linear bias set at a fixed voltage, corresponding to ap- 
proximately 1/2 maximum pulse height for a neutron energy equal to the 
mean energy of a complete run. In this example, the bias would be set 
at about pulse height 70 ^ l>ut the nan could not extend as low as 0.28- 



Mev neutron energy. 



Number of Pulses 
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experiment on li and Fe f^(E - E ) because the lab recoil 

D K 

energy of the residual nucleus was negligible. Thus the ratio 

f jE)/f JE - E^) in Eq. (29) was ®* 1 . For the lithivun and iron experi- 
D D R 

merits ve determined this ratio by comparison of pulse height analyzer 
spectra of linear pulses taken during "calibration" and "effect" runs. 
We simplified the ratio f (E, R )/f [(E - E ), R ], vhich 

1C C 1C x\ 1C 

is a ratio of intrinsic crystal efficiencies, in a reasonable way, by 
assuming 'vdien necessary (in lithim and iron experiments) that the 
energy dependence of f^ was separable from the distance dependence, 
R^. Recalling the calibration equation, Eq. ( 28 ), we see that 



C^/[Y(E) fj^(E)] = f^(E,R^) 



(30) 



Thus, if we divide the plot of versus E by the incident flux, 

Y(E), obtained from monitor neutron counts, and by fj^(E), obtained 
from linear pulse spectra, ve see that the energy dependence of f 



is determined. The same energy dependence, by assumption, applied to 

f^(E - hence by translation of the curve f along the energy 

axis by an amount E (calculated from kinematical relations \ie 

R 

obtained the ratio f(E - E ,R )/f (E, R ) insofar as the (separated) 

R X X C 

energy dependence is concerned, i.e.. 



f (E, R )/f [(E - E_), R ] S €. (E) f (R )/f (R ), (3I) 

X c'' x ' R^'’ X 4' ' x' c" x' x^’ ' 



where €|^(E) is the just determined energy dependence of the f^ ratio. 
It will be useful to write 

n /n n = eJR ^/a ), 

c' X s lx' s" 



(32) 
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where A is the scatterer area projected perpendicularly to the beam 
s 

axis and is a small correction for the degree to which the three 

solid eingles are not exactly expressible as projected area/distance 
s(iUared. 

Substituting (3I) and (32) in Eq. ( 29 ) we obtain 



f (R ) p fJE) 

da(e,E) = c^(E) /A^) 



X X 



R' 



X cr^(l - exp(- n d))’^ K C^/C^ 



( 33 ) 



The quantity e, f (R )/f (R ) = €_(R ) represents the inverse 

of the effective crystal efficiency as a fvinction of distance, R^, 

normalized to its value at R = R . If we make the reasonable assump- 

X c 

tion that this quantity is very close to vinity for R^ equal to and 

greater than R (R was 28 in. in our experiment) then neutron counts 
c c 

versus distance can be compared to inverse square law behavior (normal- 
ized at large distances) to obtain "Ehis measurement was made 

using a PoBe neutron source at several distances (Rj^) ' i’^om the 
crystal, to obtain relative counts in unit time for; each position. 

The results are shown in Fig. J. Large corrections would be necessary 
for distances less than 5 in. In the present experiments the factor 
CgfR^) was typically 1.10 to 1.20 (R^ varied from 5.5 to 9 in. de- 

pending upon the polar angle 0 of the cross section measurement). 

In view of the foregoing reductions we may write the most use- 
ful form of the expression for cross section: 
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_2 

Figo 7* deviation from (distance) efficiency response of 

a 2" dia X l" long stilbene crystal, oriented with a diameter parallel 

-2 

to the distance coordinate. Curve A: (Distance) fitted to the data 

with largest abscissa. Curve B: Smooth curve through the experimental 

points « Curve C: Calculated curve using graphical analysis and effi- 

23 

ciency corrections from Marion and Fowler. 



COUNTS/UNIT TIME RELATIVE UNITS 
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'3 - Eb )‘1 

X 0^(1 - exp(- n d))"^ K Cjj(E)/C^(E) (3^) 

where we have included a factor denoting the conversion of the 

cross section to the CM frame. For all except the experiments on iron 
and lithium, the factor ej^(E) fj^(E)/fp(E - E^^) was xmity, which 
simplified data analysis. 



Total Cross Sections 

The data analysis for total cross sections is described by 
a^(E(CM)) = €^(ND)"^ ln[(Cp - Cg)/(C^ - C^)] (35) 

where o^(E(cm)) is the total cross section as a function of mean 
center of mass neutron energy, is (as before) a symbol of the 

conversion of the energy dependence from lab to CM frame, ND is the 
number of nuclei per xmit area in the sample. In means natural log- 
arithm, C_ is the number of direct neutron coxints, C„ is the number 
i/ B 

of background coimts, C^ is the number of transmitted counts, and 
where C^, Cg, and C^ are all functions of neutron energy and are 
all normalized to a fixed accumulated proton charge. No monitor neu- 
tron normalization was applied because only count ratios are involved 
with the machine energy at one setting. 
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Error Analysis 



Differential Elastic Cross Sections 



Contributions to the uncertainty in the absolute magnitude of 
the differential elastic cross sections were present in each factor of 
Eq. (3^)^ reproduced below; 






X 0^(1 - exp(- n a d))"^ K C^/c^ 



i3h) 



Statistical errors entered the results through the factor 

C„/C . Those in C were negligible . Those which were in C„ are 
N c c N 

separately described in the figure captions for each set of experimental 
data. In general, the statistical errors were rather small (from 1 to 
3 percent); hence, they had little effect on the precision of the abso- 
lute differential cross section measurements. 

Experimental uncertainties in the factors of Eq. (3^) ai’© 
discussed below: 

1. Factor Cj^/C^. There were two experimental uncertainties 
in C^ . The detector was located at the scatterer "in beam" position 
during the calibration run. Its distance from the neutron source was 
about 28 in. The precision with which the crystal could be placed at 
the exact scatterer position was ±l/8 in. Hence, there was an ±1 
percent uncertainty from this source. Also, we include in an 

error associated with the precision to ■vrtiich the crystal was oriented 
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with respect to the incident neutrons . The aspect of the crystal rela- 
tive to the scattered neutrons^ when the crystal was in detection 
position (see Fig. l), should have been the same, because deviations 
produce differences in effective crystal efficiency. We assign an ad- 
ditional uncertainty of ±2 percent to Cj^, as a measure of the preci- 
sion with which the crystal aspects were adjusted. 

2 

2 . Factor ( R /A ) . The experimental uncertainty in R , 
the crystal to scatterer distance, was ±l/8 in. for mean R^ about 7 
in.; thus, R^ had an extreme uncertainty of ±2.5 percent. The dimen- 
sions and aspect (with respect to the beam axis) of the scatterers were 

not precisely known, which led to a ±1 percent uncertainty in A^. A 

2 

total ±6 percent uncertainty was assigned factor R /a . 

^ s 

3. Factor cr^(l - exp(- n d)) Uncertainties in the value 

of enter this factor only in second order terms in n d. The 

T T ' 

precision to which the factor nd was known determined the experimental 
error. A ±1 percent uncertainty was assigned this factor. 

k. Factor K. Althoiogh the ciorrent integrator used in deter- 
mining the normalization to fixed proton charge is rated at ±1 percent 
reproducibility, the mean normalization for several co\mting periods 
of the same accumulated charge was quite precise. So also was the mean 
normalization for several coimting periods of a different fixed proton 
charge (i.e., the counting period for calibration points was much less 
than that for riins involving detection of scattered neutrons). The 
\mcertainty in K was at most ±0.5 percent. 
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5. Factor e„(R ). This factor was obtained by comparing ex- 

-2 

perimental points with a theoretical (distance) behavior for effective 



efficiency versus distance (see Fig. 7)> Figvire 7 shows that the data 
did not precisely fix the normalization of the curve at large distances. 
The figure also shows some discrepancy between experimental and calcu- 
lated points (discussed elsewhere in this report). Considering R =7 
in. typical, we estimated ±3 percent imcertainty in 



sions involved in quoting results in terms of mean CM energies, CM 
angles, and CM intensities, where CM is center of mass. Quantities 



dominate these conversion formulas, the uncertainty in lab angles and 
energies playing a minor role. An extreme uncertainty ±0.5 percent 



7. Factor e,,(E) f^(E)/f (E - E^). In all except the lithium 
4 U D K 

and iron experiments, this factor is unity to good precision. In 
lithium and iron cross sections it is a source of additional uncertainty. 
These vincertainties, which are potentially large for a single data 
point, are reduced if one observes a smooth trend in the factor when it 
is determined for points at many energies. We estimated a ±2 percent 
uncertainty for iron cross sections and ±5 percent \oncertainty for 
lithium cross sections from this factor. 




Mj^/Mg and 




which are precise in the present usage. 



was assigned factor €^. 



Some of the experimental errors, given in 1 through 7» were 

2 

correlated; for example, R is anti-correlated to 6 o(R )• For the 

X ^ X 



purposes of this experiment, detailed examination of such correlations 
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vfas not necessary and the procedure adopted was to assume independence 
of the errors; thus, the sum of the above percent uncertainties applied. 

We concluded that exclusive of statistical errors and those 
errors associated with multiple scattering corrections (see Appendix 
E), the absolute differential cross sections of zirconium, bariimi, and 
cerivim were known to within ±l4 percent; of iron to within ±l6 percent; 
and of lithium to within ±19 percent. 

Total Cross Sections 

Measurements of total cross sections were accomplished in a 
neutron background of less than 5 percent of direct neutron counts which 
minimized second order background corrections to the data. The samples 
had small dimensions perpendicular to the incident beam direction, be- 
ing just sufficient to completely shield the detector from direct 
neutrons. The sample-to-detector and sample-to- source distances were 
large compared to the detection crystal dimensions. Neutron transmis- 
sion throiagh the samples was approximately 60 percent of the incident 
flux. All these are excellent experimental conditions for measurement 

2k 

of total cross sections. A crude calculation indicated that in- 
scattered neutron corrections to the cross section were at most 2 per- 
cent. The main constraint on the amoimt of in-scattering was the 
small solid angle the detector crystal presented to in-scattered neu- 
trons. An uncertainty of ±3 percent was assigned to accoxmt for second 
order background scattering effects and sample in-scattering corrections 
which we did not calculate in detail. Thus, we are assviming in effect 
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that the uncertainty in our rough calculation is larger than the cor- 
rection. A ±1 percent extreme error was estimated for the measurement 
of ND, the number of scattering nuclei per unit area. 

Hence, exclusive of statistical errors which are given sepa^ 
rately for each experiment, the total uncertainty in our absolute total 
cross section measurements was ±4 percent. 



CHAPTER V 



EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 

Energy and Angle Dependence of the Cross Sections 

Although the basic purpose of the experiment was to search for 
threshold effects in neutron cross sections, data such as angular dis- 
tributions and differential elastic and total cross sections were 
obtained which are interesting in their own right, exclusive of possible 
cusp effects. These data will be discussed briefly, with emphasis on 
illustrations to clarify the discussions. 

The measured angular distribution of 1.6-Mev neutrons scattered 
elastically from natiaral cerium is presented in Fig. 8. The data 
points are cross plotted from measurements as a function of energy at 

separate angles. As shown in the figure, an optical model prediction 
17 

of Moldauer, modified for multiple scattering and finite scatterer 

effects (see Appendix E) agrees very well in shape and, except at back 

angles, in magnitude, with this data. 

Figure 9 shows the measured angular distribution of 600 kev 

(cm) neutrons elastically scattered from natural lithi\mi. The data of 

this experiment compares favorably with that from a compilation of 

21 

other experiments by Howerton, which is also shown. This substan- 
tiates our estimate that multiple scattering corrections were small for 
the lithiiffli experiments (see Appendix E) . The shape of the distribution 
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Fig. 8. Angular distribution of neutrons elastically scattered 

17 

from natural cerium at = 1.6 Mev. Curve 1: A Moldauer optical 

model prediction. Curve 2: Curve 1 modified with a multiple scattering 

correction. Curve 3t Curve 2 modified with calculated experimental 
finite scatterer effects. Triangles are experimental points from the 
present experiment . 
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Fig. 9“ Angular distribution of neutrons elastically scattered 

from natural lithium at E = 0.6 Mev. Open circles are data from the 

n 

present experiment. This data is not corrected for multiple scattering. 

27 

Triangles are data from Howerton's compilation. 
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is interesting, being peaked backwards, presxjmably as a result of inter- 

28 

ference between the strong £ wave resonance at 0.25 Mev and an s 
wave background. This backward peaking means that one should look at 
forward angles for cusp effects, because the theoretical percentage ef- 

1/2 

feet is inversely proportional to (da) ' (see Chapter II). 

In Fig. 10 ve show some of the differential cross section data, 
versus energy, from which the angular distribution. Fig. 9 , was cross- 

28 

plotted. At the lower energies the upper edge of the £ wave resonance 
at 0.25 Mev is in evidence. The wavy structure in all the cross sec- 
tions is probably not an experimental effect; it is q\xite surprising, 
because for a light nucleus such as lithium one might expect rather 
smoothly varying cross sections at these energies, except for isolated 
resonances. The inelastic cross section to the first excited state of 

7 29 

II , as measured by Freeman, lane, sind Rose shows similar structure. 

The differential elastic cross section of iron, versus energy, 
was measured for four angles using a rather thick neutron target 
Table l) in the vicinity of 0.8^7 Mev. The res\ilts are shown in Fig. 

11. It was later determined that a much thinner neutron target was 
needed to resolve the struct\ire and to search for cusps (see Fig. l^t-; 
target Table I was used). 

Energy and Angle Dependence of Threshold Features 

We discuss here the experimental results as they pertain to 
Wigner and "energy-averaged" Wigner cusps. Frequent reference will be 
made, explicitly and implicitly, to the predictions of cusp theory. 
Chapter II. 
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Fig« 10. Neutron differential elastic cross sections (CM) of 
natural lithium versus center of mass (CM) energy for five angles (CM)o 
T he data are not corrected for multiple scattering. The original 
experimental results have been averaged such that statistical errors 
in these points are approximately ±1^. 
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Fig. 11. Neutron differential elastic cross sections of natural 
iron versus energy for four angles. The data is not corrected for 
multiple scattering. Incident neutron energy spread is approximately 
25 kev. Statistical errors are indicated by the size of the symbol. 
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Lithium (92.6^ Li^) 

The cusp results for lithium were essentially negative. The 

theory of Appendix B, Case I, applies. We saw in the preceding section 

that the largest percentage cusp effect in da is expected at forward 

angles. Figure 12(b) shows that an experimental effect was indeed 

present at (CM) in the vicinity of O.L 78 Mev, the first excited 

7 

state of LL . However, using the reaction cross section magnitude of 

29 

Freeman, Lane, and Rose, and the applicable theory, we find that the 
observed effect was at least a factor 5 too large to be explained as 
a Wigner cusp; an effect of approximately 2 percent is predicted 50 kev 
from threshold. Figure 12(a) compares our measurement of lithium total 
cross section with that from Howerton's compilation. Although the 
two sets of data agree, no Wigner cusp is in evidence in either. We 
had hoped to see an effect in the total cross section, because the re- 
sults of Lane, ^ (see Fig. I3) appear to show a break near O .5 

Mev in the coefficient of the cross section expansion (see Chapter 
II). The energy scale change in Fig. I 3 at about O .5 Mev may bias this 
interpretation of their data. 

Iron ( 91 o 7^ Fe^^) 

With iron distinct features appear in all cross sections at the 
energy of interest (2 first excited state Fe^ , O.8U7 Mev). This is 
exemplified in Fig. iL where a resonance feature appears, apparently 
at an identical energy (just above 0.847 Mev) in the total cross sec- 
tion, in the l45“ differential elastic cross section and in the reaction 
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Fig. 12. Neutron cross sections of natural lithium versus 
energy, (a) Total neutron cross section; triangles are data from 

27 

Howerton's compilation and solid circles are the energy-averaged 
(a E =20 kev) data from the present experiment, (h) CM differ- 
ential elastic cross section of the present experiment. 
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by Li"^. 



Fig. 13. Differential cross section for scattering of neutrons 

Reproduced from the work of A. Langsdorf, Jr., R. 0 . Lane, and 
20 



J. E. Monahan. 
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Fig. 2. The differential cross section for scattering of neutrons by Li^. 



COS 6 



77 



Fig. l4. (a) Neutron differential elastic (l45“) cross section 

of natural iron versus energy, (b) Total neutron cross section of 
natural ironu (c) Reaction cross section Fe^ (n, n*, x (-847 Mev)) Fe^ . 
Statistical errors are the size of the symbols in (a) and (b), but are 
±10^ in (c)o The smooth curves show the result of removing a calculated 
Wigner cusp effect londer one possible assumption on the phase of the 



scattering amplitude. 
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cross section Fe^^(n, n')Fe^^ (O.8U7 Mev) . The latter cross section 

19 

was measured by Tucker, et using the same tritium target. An- 

other resonance feature which seems to be a "mirror image" of the 
aforementioned resonance across threshold is evident in these data, 
and in all the differential cross sections at other angles investi- 
gated [see Fig. 15 (six angles)]. Application of the most general 
theory is hopeless because of the large number of unknown phases. At 
least six phases contribute to the cusp, even with a restriction on 
participating incident partial waves to those with orbital angular 
momentiim i less than 3 ; see Appendix B. A reasonable and tractable 
special case of the theory, which might apply to the iron data is dis- 
cussed and explicitly evalviated in Appendix B, Case II. Unfortunately, 
the special assumptions involved allow cusp analysis only over a nar- 
row energy region, entirely above threshold. The resiilts of this 
analysis are superimposed on the data in Fig. 15 for differential cross 
sections (6 angles) and in Fig. 16 for the total cross section. The 
results are neither startling nor conclusive, although by removal of 
the cusp effect the peak energies and widths of the above threshold 
resonance in the differential elastic, total and inelastic cross sec- 
tions can be broiaght into exact coincidence. 

Ceri •urn- Bari um-Zirconim 

We consider now the results obtained for targets of higher 
atomic number in the region where the optical model is expected to ap- 
ply because compound resonances are closely spaced — cerivmi, barium, 
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Fig. 15. Neutron differential elastic cross sections of 
natural iron versus energy for six angles. The data is not corrected 
for multiple scattering, (a) gives the data over the full energy 
range of the rim. The lined box shows the section of the data which 
is expanded into (b). In (b) statistical errors are indicated by 
the symbol size. Open triangles show the result of removing a cal- 
culated Wigner cusp effect, under one possible assumption on the 
phase of the scattering amplitude . 
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Fig. l6. Total cross section of natural iron. The size of 
statistical errors is indicated by the size of the symbol. Triangles 
show the result of removing a calculated Wigner cusp effect \mder one 
possible assumption on the phase of the scattering amplitude. 
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and zirconium. The theory of "energy-averaged" Wigner cusps is then 
also expected to apply (see Chapter II). 

Composite data for cerium ( 88 . 5 ^ Ce^^), given in Fig. Yf , is 
representative of the cusp analyses we present for all three elements . 

A special feature of this figure, however, is that we present the 
original data, to show the degree of scatter on interleaving mins (most 
of which is not statistical, but experimentally reproducible). In the 
remaining figures presented for this group, the data was averaged over 
an energy interval suitable to reduce statistical errors to about ±1 
percent (refer to Chapter IV). This procedure standardized data presen- 
tation for a large number of angles. 

This and other figures in this group also show cross section 
trends predicted by the optical model. These were obtained from the 

30 

work of Campbell, et The optical model shapes were interpolated 

:^om very few points given in the reference and should not be consid- 
ered a precise computation; hence, they were nomalized to the experi- 
mental data at threshold energy. The magnitudes of the predicted cross 

30 

sections were consistent with the estimated experimental errors 
(Chapter IV) and the corrections given in Appendix E. 

A cusp effect is quite apparent in Fig. 17(b) by comparison of 
the optical model trends with the experimental cross section. To com- 
pare the magnitude of the cusp effect with theory one could subtract 
the experimental cross section from the optical model prediction and 
plot the difference as a fiinction of neutron energy. But, as mentioned 
before, the optical model trends are somewhat uncertain and, furthermore. 
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Fig. 17- The experimental inputs to a ceriiom "energy-averaged" 

Wigner cusp analysis, (a) Energy-averaged (A E = i»0 kev) experimental 

ll(Q 2.9 

inelastic cross section Ce (n, n' y (1.597 Mev)) Ce . (b) Origi- 

nal data for elastic 95" scattering of neutrons from natural cerium. 

The two symbols indicate the reproducibility of the data on successive 
interleaving energy runs. The size of the symbols indicates the size 
of the statistical errors, (c) Same data as (b) corrected, point by 
point, for a calculated "energy-averaged" Wigner cusp. The smooth 
curves in (b) and (c) are optical model predictions of Campbell, et al. , 
normalized at 1.6 Mev. 
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the experimental cross section has natural fluctuations, which will he 
discussed below. Hence, it was thought more appropriate to remove from 
the experimental cross sections the theoretical cusp effect and to show 
(a) that in the resultant curve there would be no more cusp effect--at 
least within the natural fluctuations of the cross section, and (b) 
that the resultant curve would follow the optical model trend, which 
theoretically represents the cross section without threshold effects. 

The procediire for removing a theoretical cusp from the measured 
cross section was to mxxltiply the reaction cross section (i.e., inelas- 
tic neutron scattering cross section) by the isotopic percentage of the 
isotope of interest, by the appropriate multiple scattering and finite 
scatterer correction (the reciprocal of Table III (Appendix E) values), 
and by the theoretical angle dependent factor (Chapter II) . This neg- 
ative quantity was subtracted from the experimental data [see, for 
example, Eq. (26)]. 

Cerium (88.9?^ Ce^^) 

We now consider the specific case of cerium, for which Eq. (26) 
applies. Tucker, et ^., measured the reaction cross section [Fig. 
17(a)] Ce^^(n, n')Ce^^ (I.60 Mev).^^ As previously mentioned, inci- 
dent d waves dominate the reaction to at least I50 kev above threshold, 
hence i' = 2 in Eq. (26). The results of removing the theoretical 
cusp effect are shown in Fig. 17 (c). There is little doubt that the 
theoretical "energy-averaged" Wigner cusp, with i' = 2 Is a satis- 
factory description of the change at threshold. Note Eq. (26) predicts 
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that if incident £ vaves were imagined to dominate this reaction, rather 
than the actual d waves, the resulting expected cusp, at polar angle 
90“, would be a factor 20 /l 3 or about 50 percent larger. The plot cor- 
responding to Fig. 17(c), corrected for an expected = O) cusp, 
would show an upward break at threshold and would follow neither the 
below threshold data trend nor the optical model trend. Thus we can 
conclude the reaction is not dominated by incident £ waves. This type 
of information can be useful in assigning excited state spins, for 
example . 

The differential elastic cross section of ceriim was measured 
at ten angles, versus energy. These results are shown in Fig. I8. 

Thej points obtained after removing a calculated cusp are also given 
for each angle. In all cases, the adjusted data so obtained follow 
the trend of the below threshold data and the fluctuations in these 
is the same above and below threshold. In general, the optical model 
curves agree nicely with the modified cross sections; that is, they 
follow what the theory defines as "cross section without reaction." 
There is no reason to expect this agreement to be perfect, since the 
optical model was constructed to explain gross structuires only, over 
a wide range of atomic number and energy. We should recall the assump- 
tion, made for Eq,. (26), that the fraction (called R^i in Chapter II) 
of the reaction probability for proceeding from compound nucleus states 
of a particular J' is the same as the fraction of the dominant inci- 
dent partial wave states (all with the same orbital angular momentum 
£') which have J = j'. Other assimiptions for Rji 



lead to cusp 
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Fig. l8. The neutron differential elastic cross section of 
natural cerium versus energy for ten angles. The data is not corrected 
for miiltiple scattering. Statistical errors are ±1^. The points are 
an energy average (A E =50 kev) of the original experimental data. 

For each angle the left side is the measured cross section. The right 
side is the cross section with a calculated "energy-averaged" Wigner 

cusp removed. The calculation is based on the excitation of the 1.6 

l4o + 

Mev state in Ce . No cusp contributions from the 0 (l.90 Mev) or 

(2,08 Mev) states are included. The smooth curves are optical . 

model predictions of Campbell, et normalized at 1.6 Mev,^*^ 
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magnitude predictions differing from those used here by as much as 30 
percent, at certain angles. Although not the case in the present work, 
it is possible that a more refined experiment could be used to measure 
R ,j such a measurement could assist in the understanding of neutron 
nucleus interaction mechanisms, particularly if it turned out that a 
particular reaction preferred one value of compound nucleus J' to 
another, even though resonances were closely spaced and possibly over- 
lapping . 

Barium (71.7^ Ba^^^) 

The measTored cross sections of barium and the results of the 

appropriate cusp analysis are compiled in Fig. 19. The reaction cross 

section Ba ^ (n, n')Ba (1.427 Mev) was measured by Tucker, et al . 

Again the spin of the first excited 1.427 Mev level is 2^ and the ground 

state 0^. That £ wave neutrons dominate the reaction very far from 

threshold is not so clear in this case, since the reaction shape is 

distorted by resonance structure. However, optical model calculations 

l8 

show that barium is near a meiximum in the _s wave strength fimction. 

As in cerium then, the reaction is considered to proceed by the Z' =2 
incident partial wave. The results of theoretical comparison are sat- 
isfactory. The cusp effect is smaller than in ceriijm because of the 
smaller isotopic percentage and the smaller reaction cross section near 
threshold. In Fig. 19(c) we see that it would be difficult to distin- 
guish a cusp at all were it not for a comparison with the optical model 
trend \diich happens to be rising. 
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Fi^;* 19- Neutron differential elastic (90° and 125°) cross 

'cc^ion,. of iiatural bari 'm versus energ^^. (a) Energy-averaged (A E = 

20 l;c.) inelastic cross section for Ba^^^(n, n' r (1.427 Mev) ) 

Gtatistlcal error.; arc ^yfo. Tlie left sides of (b) and (c>) are an 

cr.cr^, a err,^ e (A E - 2y !:e /) of the original expex-imental data. Tnis 

datr. is :ot correcocd for Multiple scattering. Statistical errors are 

±l/j. The right sides of (b) and (c) shov the results of rer.ioving 

-al^-ulcAcd "en.crgy.-av’eraged" Wigner cusps. Tlxe smooth cur\es arc op- 

, 30 

deal ; wdel predictions of Campbell, et al., normalized at 1.427 Mev. 
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ZTrconium 

The measured cross sections of zirconium and the results of 

cusp analysis are compiled in Fig. 20. The inelastic reaction 

Zr^^(n, n')Zr^^*(l.75 Mev) vasj measured by Tucker, et The spin 

of the 1.75-Mev first excited state is 0 and the ground state 0 • 

Tucker also showed that £ wave outgoing neutrons dominate the! reaction 

not far above threshold. Hence we assume that the reaction proceeds 

90 

from incident p waves, i.e., i* =1. For Zr^ the low isotopic per- 
centage and small reaction cross section near threshold combine to 
produce a small calculated cusp effect. At 60'" polar angle the pre- 
dicted cusp effect is too small to be distinguished in the data. Note 
also that at 6o® the optical model trend does not agree well with the 
measured cross section slope. At 130 the cusp is distinguishable and 
the theoretical cusp fit is satisfactory « 

Total Cross Sections 

The theoretical “energy-averaged" Wigner cusps which were used 

to explain the preceding experimental results (Ce, Ba, Zr) are based 

22 

on the same assumptions that underlie the optical model <> These as- 

s-umptions also lead to predictions that the energy average of the total 
cross sections should be smooth across thresholds [see Eq. (23)]. Our 
measurements of the total cross sections of cerium and barium are given 
in Fig. 21. The predicted cusp effect for the total cross section of 
zirconium — even if the energy-averaged theoiy would not apply — would 
be so small that we did not investigate this case. For completeness. 
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Fig. 20. Neutron differential elastic (6o® and 130°) cross 

sections of natural zirconi\jm versus energy, (a) is a ^40 kev average 

90 90 19 

of the reaction cross section Zr (n, n') Zr (1.75 Mev) . Statis- 
tical errors are ±55^. The left sides of (b) and (c) are an energy 
average (A E =50 kev) of the measured cross sections. This data is 
not corrected for multiple scattering. Statistical errors are ±1^. 

The right sides of (b) and (c) show the results of removing calciilated 

"energy-averaged" Wigner cusps. The smooth ciirves are optical model 

30 

predictions of Campbell, ^ aX., normalized at 1.75 Mev. 
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Fig. 21. Neutron total cross sections for natural bariijm, (a); 

cerium, (b); and zirconium, (c). Statistical errors are indicated by 

the size of the symbols. The data on barium and cerium were measured 

in the present experiment. The data on zirconium is taken from 

27 

Howerton's compilation. In (b) circles are the total cross section 
of cerium; triangles are integrated elastic cross sections, corrected 
for multiple scattering; and squares are the s\jm of the integrated 
elastic and of the reaction cross section (corrected for isotopic per- 
centeige) Ce^^(n, n' y(l.6 Mev))Ce^^ » 
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27 

the zirconium total cross section data from Howerton is also presented 
in Fig. 21, showing the inadequacy of widely spaced energy points for 
the present study. The data for barium and cerium were taken with an 
energy spread of about 30 kev, yet they are certainly not smooth f\mc- 
tions of energy, displaying marked structure » The mean values of these 
fluctuations do not seem to drop markedly as threshold is crossed, sug- 
gesting agreement with the theory o If, however, the theory of energy- 
averaged cusps does not apply, but one assumes instead a pure Wigner 
cusp phenomenon, the cusp predicted is at most a 5 percent effect which 
would be masked by the fluctuations. Hence, the total cross section 
data cannot be considered a verification of either assumption. 

That the fluctuations are not experimental is supported by the 
differential cross section data for cerium at forward angles (Fig. l8) 
which shows the same structure . Note that the structure does not con- 
tain recognizable resonance shapes. It is possible that the consider- 
ations of Ericson on statistical fluctuations in cross sections might 
31 

apply here. He shows that if a compound system is characterized by 

many overlapping resonances, each of width P, and with relative ran- 

21 

dom phase of the value quantities y of resonance theory, cross 

K 

sections would have a fluctuating structure of average width, r. 

In the case of cerium only, sufficient differential elastic 
data was taken to permit integration over all solid angle. The total 
elastic cross sections so obtained for several energies, after correc- 
tion for multiple scattering, are plotted in Fig. 21(b). The drop at 
threshold (1.6 Mev) is evident, as, in accord with any theory of 
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threshold effects, flux conservation alone requires the total elastic 

cross section to drop above threshold for a reaction. In the case of 

lUO 

Ce , the experimental differential cross sections all drop, so the 
total elastic cross section does also. It is reass\iring that, when 
the inelastic scattering cross section is corrected for isotopic per- 
centage and added to the integrated elastic cross section, the result- 
ant cross section (ioe., the total cross section) is smooth across 
thresholds The agreement in magnitude between the latter curve and 
the directly measured total cross section is considered somewhat 
fortuitous . 



CHAPTER VI 



CONCLUSIONS 

The experiments have shown that under suitable conditions the 
threshold effects expected in neutron elastic scattering can be ob- 
served. From an experimental point of view particular attention must 
be paid to the choice of energy resolution consistent with observing 
a cusp effect in reasonable time . 

Explicit theoretical expressions have been obtained for ex- 
pected Wigner and *’energy-ave raged*’ Wigner cusp effects under general 
angular momentum conditions^ within some reasonable restrictions imposed 
in the derivation. The effects were found to be contained in the coef- 
ficients^ of the Blatt-Biedenharn formulation of the general 

scattering problem. Since many survey experiments involving neutron 
scattering have been analyzed in teims of the B^^ coefficients^ we 
make the obvious siiggestion that such results be examined near excited 
state energies for possible tests of cusp theory. 

The expected Wigner cusp in lithium cross sections was expected 
to be an approximately 2 percent effect^ 50 kev from threshold. Our 
experiment was not sufficiently refined to pick out an effect of this 
size . 

Resonance structiire in the measured cross sections of iron made 
a Wigner cusp interpretation ambiguous. Features occurred at the thresh- 
old energy^ but their analysis in terms of cusp theory was inconclusive. 
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The measured differential elastic cross sections of Zr, Ba, 
and Ce dropped at inelastic thresholds. Analysis of the drops in 
terms of "energy-averaged" Wigner cusp theory resulted in good agree- 
ment between theory and experiment. This indicates that the usual 

energy-averaging assvimptions underlying the optical model are indeed 

7 

applicable to these nuclei . Also, the cusp theory shows that there 
cannot be a large amount of direct interaction between incident neu- 
trons and these nuclei at the neutron energies used in these experi- 
ments 2 Mev), or else below threshold effects would also be seen. 



APPENDIX A 



THE PHYSICAL BASIS FOR THE S-MATRIX AND ITS USE 
IN CROSS SECTION FORMULAS 

We will review the principles of scattering analysis and the 

12 

definition of the S-matrix as presented by Blatt and Biedenharn. 

They pointed out that general expressions for cross sections in terms 
of the S-matrix are rarely given in standard texts . 

The general form\ilas are applicable to any collision process 
in which two particles collide and two particles emerge. However, in 
this appendix, and after a certain stage in the development, only the 
special case of elastic scattering with no i change will be con- 
sidered . 

The channel index a defines the type of incoming particle and 
the state of the struck nucleus. The channel spin, s, is the total 
spin angular momentum of the channel, obtained by vector addition of 
the intrinsic spin, i, of the incident particle and the spin, I, 
of the struck nucleus. The orbital angular momentum (in the center of 
mass system), Z, is combined with s to form the total angular 
momentum, J. Only in this particular appendix, primes will be used 
to refer to outgoing or final state quantities. J is conserved during 
the collision, hence the values of Z' and s' must be consistent 
with J' = J. Similarly, the parity it, of the initial system, is 
conserved, and it = it ' . 
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The probability amplitude Si f.i „ for a collision with 

CK S i/ j CtS Jo 

total angular momentum J, from channel Q!si to channel Ot's'i' is 
referred to as an element of the scattering matrix. The scattering 
matrix, or S-matrix, is defined by considering asymptotic forms of the 
collision wave functions, as follows. 

Consider one definite value of J and = M. At any particu- 
lar energy, E, the channel wave number, k^, and relative speed, 
v^, in each channel a are given by the energetics of the reaction. 
Only those channels which are "open" (k^ real) are considered. Let 
9 ^ be the product wave f\mction of the target nucleus and incident 
particle. The spin and angle dependence of a wave function of total 
angular momentimi J, = M, orbital angular momentimi i, and chan- 

nel spin s is given by 




i s 

= E S ( -^sm 

m.=-i m =-s 
£ s 







X 

s,m 
’ s 



(Al) 



where X is the incident spin wave function and Y. is the well 

s,m ^ £m 

s 

known spherical harmonic of angular momentum £ and projection m. 
Using these definitions, the most general wave function in channel CC, s 
with total angular momentum quantum numbers J, M, consists of the 
superposition of incoming and outgoing spherical waves. Asymptotically 
(at sufficiently large distances) we can write (for neutral particles) 






_ JM r . / . 

B . expL+i(k r - 

as^ ^ ^ a a 



I £ n)]| (A2) 



Ill 



Substitution of Eq. (A2) into the wave equation^ yields a unique connec- 
tion between the outgoing wave coefficients, B, with the ingoing wave 
coefficients, A. Thus 



JM 



B„_ I „ I I — 2 2 2 I _ I 0 1 

OC s £ 



J „ JM 



a s £ a s i ,*asi ocs£ 



(A3) 



Eq. (A3) defines the scattering matrix, With these definitions it 

is evident that ^ is unitary (conserving the wave function normaliza- 
tion in time), symmetric (reciprocity theorem), and independent of M. 

It remains to decompose the incident plane wave into the f mo- 
tional forms we are using. The incident wave is cp exp(ikz)X 

cts s • m 

^ s 

First, one writes 



m 



where 



i+s J 

Y, ^ = 2 2 (isom^lis 



M 



J=U-s| M=-J 



Jis 



Using the asymptotic form of J’». i for kr > > i the following ex- 

2 

pression is obtained for the plane wave in channel a, s with spin 

direction given by m ; 

s 



- / \f- w J J+S 

exp(ik z )X CD ^ ISJLL cp E Z Z 



a a' s,m ^as k r^ ^as 

^ c’ r>t 



a a J=0 M=-J ^=|j-s| 

X (^gsom^jis JM) i^{2£ + l)^ ^ 



M 



J£s 



(aU) 



X |^exp[-i(k^r^ - ^ £ n)] - exp[+i(k^^ - i ^ «)]|- 
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Comparing Eq. (A4) vlth the standard form, Eq. vre see that the 

amplitudes of the ingoing spherical waves are given by 



^asf = (isomjis JM) i^(2i + 1)2 



(A5) 



The amplitudes of the outgoing spherical waves are determined from Eq. 
(A3). As in the main text, we now specialize to elastic scattering 
with no change of £ or s. They are 



JM 



1 J+s 



J-sl 



X (fecm^l^s JM) l\21 + 1)= 



(a6) 



In order to investigate cross sections ve only need that part of the 
plane wave solution^ Eq, (a 4)^ which is due to the scattering interac- 
tion: i.e,, we want to obtain ^ . from \|r = _ . + 

^scattered ^ incident 

\|r , Therefore, we substitute Eqs. (A5) and (a 6) into Eqo (A2) 

scaTr uereci 

and subtract the incident wave, Eq. (a 4); the result is 



1 



00 J J+s 
ZEE 



(A7) 



, , . = 1 ^ — cp 

^scattered a r^ ^as », t « i t i 

a J =0 M=-J i=|J-s| 

P i 

X (isom I is JM) i (2i + l)2 
s 

X exp[l(k^^ - i i „)1 X [1 - 

If the spin angle function is decomposed using Eq. (A1), then 

Eq. (A7) can he used to define the scattered amplitude, 
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rjis 
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■where q is explicitly 



00 j J+s I 

M-J i=lj-s| K=-i 



.1 i 

X TT^(2i + l)^(isom jis JM)(isMm' |is JM) 






(A9) 



The differential cross section for the scattering asm^ 
given hy 

O p 

da , = ^ q„ I 

asm : asm a ^sm ;asm 

c c c C i 



asm 1 s 
s 



(AlO) 



For unpolarized incident beams, we must average Eq. (AlO) over initial 



spin direction m and svan over final spin direction m ' : 
s s 



da 



OSjOS 



= (2s + 1) 



-1 



s 

Z 



s 

Z 



da 



m =-s m =-s 
s s 



asm' : asm 
s*^ s 



(All) 



This expression is greatly simplified by sum rules for Clebsch-Gordan 
coefficients, and it is from this calculation, as performed by Blatt 
and Biedenharn, that one obtains the useful final cross section expres- 
sions, involving Z coefficients. 

Finally, to obtain the differential elastic cross section, one 
must average over the possible values of incident channel spin, s. 



, - 2s + 1 

a,*a s=|l-i| (21 + l)(2i + 1) ^'^as;as 



(A12) 



Equations (A12) and (All), as explicitly vorked out in Ref. 12, are 
just Eqs. (l), (2) and (3) of the main text. 



APPENDIX B 



SPECIAL CASES OF INTEREST TO THE PRESENT EXPERIMENT 

Case I. Cusp effects in neutron scattering from Li 
near the threshold for excitation of the first 

^ n 

excited (O.U 78 Mev) l/2 state of LL (grovind 
state spin and parity, 3/2 ) • 

The large resonance in the elastic cross section at about 0.25- 

Mev neutron energy has been analyzed and found to correspond to a 3 

8 28 *7 *7"^ 

state in Li . The reaction cross section Li (n, n*)Li ( 0.478 Mev) 

has a broad resonance at 1 . 35 -Mev neutron energy, probably associated 

+ 29 

with a compound 1 state. Also, the inelastic scattering is domi- 
nated near threshold by outgoing £ waves, with compound angular momen- 
Tt - 29 

turn J = 1 . Finally, any process involving incident partial waves 
with ^ > 1 is unlikely at 0o5 Mev because kR 0o3* We conclude 
that in the vicinity of 0.48-Mev incident (center of mass) neutron 
energy only s and £ waves participate significantly in elastic scatter- 
ing. Also, incident s waves dominate the inelastic scattering. The 

latter conclusion results from the known spins and parities of the 

7 

ground and first excited state of II , and the assumption that no 
process with £ > 1 is likely o With s waves dominant in the outgoing 
reaction channel the only possible incident channel leading to the re- 
action is characterized by J* = 1, s’ = 1, and £' = Oo Using these 
values, and I = 3 / 2 , i = l/2, in Eq. (12), one obtains 
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A B (1;1) = - J s J 0 1, 1 L) 

^ J£ 



X 2 sin 6 
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1£ 




sin(26Q^ 

cos(26q^ 




(Bl) 



Restricting ourselves to £< 2, as indicated above, we need these 
Z coefficients 



Z^(0101,10) = 3 , z^(ll01,ll) = 3 , Z^(1201,11) = 5 



2 

and all other Z (iJOl^lL) values are zero. Therefore, 

2 sin^ 

A Bq( 1;1) = - I 1 • pel } 

'^sin 2 &q^ 

sin(2S^, - 6^ ) 



AB^(1;1) = -(3/2)a^^ 



sin 6 



5 2 
1 sin 6^^ 
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:os(25ji - 
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(B2) 



Equations (B2) might be compared to data such as that of Langsdorf, 

20 

Lane and Monahan, reproduced herein as Fig. I 3 . Unfortunately, the 
energy scale in this figure changes in the region of interest ( 0.5 Mev 
lab)o We present the figure primarily to illustrate the utility of 
cusp expressions in the fom of Eqs. (B2). The cusp prediction for 
the differential elastic cross section is obtained by substituting Eq. 
(B2) in Eqs. (9-H) to obtain 
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(B3) 



Using Eq. ( 16 ), we obtain the total cross section cusp term, 

,1 

^ (B4) 

"01 



r cos 26 'I 

Ao,-(3/8)(c^,^) I “} 

- Sin 25^ ~ ^ 



vhere (a . in Eqs. (B3) and (B^) is the experimental reaction 

^ ri =0 exp 

cross section near threshold. 

7 

For the case of LL ve recall the restriction imposed on the 
derivation leading to Eq. (12), that the amount of scattering with 
change of s, hut no change in energy or orbital angular momentum Zy 
be small. There are no good reasons for believing this to be true in 
the present case. Change of ? = ? + i^ from |s] = 2 to |s^| = 1 
(s^ means outgoing "elastic" channel spin), or vice versa, in the li 
case is essentially a spin flip followed by a recoupling of s and 
i( <6 g l), and this is just the type of interaction involved in the 
reaction leading to the first excited state. We have seen, however, 
that in fact £ waves {&' = O) dominate the LL reaction; thus "spin flip" 
elastic scattering (which requires i 5 l) does not couple to the reac- 
tion. Hence, Eqs. (B3) and (b 4) apply. Using the reaction cross section 
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29 

value ( 100 kev above threshold) of Freeman, Lane, and Rose, and the 
present experimental values of do (3^° CM), we can estimate the ex- 
pected percent effect to be 

A do/da* = (3/32 Jt) a^/da® s (3/32 «) .05/.05 s 3 percent effect 

(~ 100 kev from threshold) 

Correspondingly, for the expected cusp in we obtain 

A = (3/8)( .O5/I0I) = 2 percent effect 

(100 kev from threshold) 

Case II. Cusp effects in neutron scattering from 
iron near the threshold for excitation of the 

4- 56 

first excited state (2 , 0.847 Mev) of Fe''^ 

(ground state spin and parity 0"^). 

Equations ( 9 ) to (12) for pure Wigner cusp effects apply to 
this case because the resonance structure in the cross sections is ap- 
parently resolved. We must compare Eq. ( 9 ), the expression for change 
in cross section, with experiment. There are too many possible phases 
involved in a general treatment of Eq. ( 9 ) to make cusp interpretation 
meaningful; therefore, we sought a tractable and reasonable simplifica- 
tion. 

The experimental data, see Fig. l4, display a resonance at an 
identical energy in the reaction cross section, in the total cross sec- 
tion, and at six angles in the differential elastic cross section. On 
this evidence, we assxmied that in the vicinity of this resonance the 
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scattering and reaction proceed through the same compound state of total 
momentum, J', and are dominated by one incident partial wave, . 

With this assumption, and substituting I = 0 , i = l/ 2 , s' = l /2 in 
Eq. ( 12 ) we obtain 



2 sin 



A Bj^(a a i) = 



Z^(i’ 






J' £' J'; 



jL 

2 



L) 



I I 

^ sin 2 6 ^ 

X, p 



T 

OVa'a(B5) 



where J ' must be either £’ + -^ or £' - depending on the nature 
of the compound resonance. By substituting Eq. (B 5 ) in Eqs . ( 9 ) to (ll), 
one obtains 



A da 



a; a 



/ J' s 
^r£ '-g- exp 

5 it 



o • 2 ^ J' 
2 sin 6 a 
^ 2 

• orJ' 
sin 26., 1 

Ji 2 



X Z i prices 6)d« (b6 ) 

Li=0 

J * 

vhere (cr .,i) is the measured reaction cross section. 
rx 2 

The optical model indicates a maximum in the £ vave strength 
function at A « 5 ^^^^ and the analysis of Tucker^ ^ shovs 

that £ vave dominance of the outgoing neutrons is consistent vith the 

56 

measured inelastic cross section in Fe . Knowing that the excited 
state spin is 2 ^^ we conclude i* =2 and J* is either 5/2 or 
3/2. For J* = 3 / 2 ; 



3/2 X 2 sin 8^' 



A da 



a; a 



(n ') 

''°r^'=2''exp 




i Sin } 






[1 + Pg (cos e)]d« (B 7 a) 
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For J' = 5/2; 



A do 



( 5/2 ) 2 sin^ 

^ ri'^^exp ^ ^ 



a; a 



i+jt 



sin 26, 



f 



8 6 ' 

1 + - P^Ccos 0) + - Pj^(cos e) 

(BTb) 



Using Eq. (l6), ve obtain the expected cusp effects in the total cross 
section. For J' = 3/2 J 

cos 25^/^ ^ 

'' (B8a) 



' <”rr^)exp 1 . 3^3/2 

22’ 



For J' =5/2; 



A CT^ = 



cos 26 






rr^'axpl. 3^^3g5/2 



^2 



} 



(B8b) 



Since ve have assmed isolated resonance behavior in the cross sections 

in this case, ve can vrite the phases 6 as explicit functions of energy » 

Defining e, = (E - E )/(r,/2), vhere E is the resonance energy 
J Cx O J O 

and r^-i is the full width at half maximum, we find 
J 

2 sin^ 6*^ = ^( 1 + e^) , 

cos 2 6'^' = - (1 - e^)/(l + e^), (B9) 

T * P 

- sin 26 =2e/(l + e ) , 

where e ranges from - » to + «> across the resonance. 

Below-threshold effects are possible in both the differential 
elastic and total cross sections if the special assumptions are valid. 
Unfortunately, this cannot be checked in the present experiment, since 
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another resonance feature appears just belov threshold, whose spin is 
unknown . 

The above-threshold feature can be analyzed for possible cusp 
effects within the framework of these assumptions, by using Eqs . (B7) 
with Eqs. (B9) substituted. 



A da 



^^ri'=2^exp ( 2 



ap. 






1 + e 



(-(1 + PgCcos 0), if J' = 3/2 



(BIO) 

dfl 



(1 + ( 8 / 7 ) P 2 (cos 9 ) + ( 6 / 7 ) Pj^(cos 0)), if J' = 9/2 

(Bll) 



A ! ^ (e - l) 

(°r^- = 2 ) exp ^77"^ 



Equations (BIO) and (Bll) shov an interesting property. The 
cusp term in expressions for elastic cross sections (differential and 
integrated) is negative over the entire resonance. The total cross 
section expression is negative for - 1 < e < 1, and positive for 
€ < - 1 and ' e > 1. This can be interpreted simply as an increase 
in the total cross section resonance vidth due to the reaction o Hov- 
ever, the experimental neutron energy spread (« 7 kev) vas too broad 
to allov measurement of such a difference in the vidths of the elastic 
and total cross section resonance featirres. 

In conclusion, ve remark that the belov-threshold feature (see 
Fig. l4, for example) vhich ve have described as "mirrored" across 
threshold is just the kind of cusp effect one might expect from the 
theory. We vere disappointed that our cusp analysis did not demonstrate 
this feature to be a cusp anomaly. 
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A PROOF 

We wish to show that for the case of many resonances in the 
averaging interval one finds under certain conditions 



sin(25'{, - 6*^ 



( a_„, 2 sin ,, t r / 

^ nn<=if PF» . - PS J 



rl' 






(Cl) 



ri" b;;, b 
0 



K _ K 
ii* ^JJ* 



,K 



where B is the Kronecker delta. Expanding the left side of (Cl), 
we obtain 



- { 



gj J 


T * 

p sin(26^, 






^ cos(26^, 


-8^)i 




J 


_ J 




Sin 62 - cos 26. 

X/ Jo 




. J 


J 


2B_g, 


Sin 26^ + sin 26^ 



) 



• 2 . J 

;in B^ 

• 2 cJ 
;in 



} 



(C2) 



For simplicity, let 



J' 


sin 26 ^, 


T, J' 1 


- COS 26 


" "ri’- 


^ COS 26^1 ^ 


= ^ “rrl 


(■ sin 26^1 



J' 

Z' 



a = sin 26 



b s 



o • 2 . J 
2 sin 6^ 
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Then expression (C2) becomes 



= ( Aa + Bb ) 



(C3) 



We now assume, for I i Z' and/or J / J' and within an averaging 
interval, that the quantities A, B, a, and b may be written as the 
sum of a mean value and a fluctuating part; for example, A = A“ + A A, 
where ( A A ) = 0 . Thus, expression (C3) may be written 

( A°a° + B^b® + + a^A A + B“A b + b°A B + AAAa+ABAb> (C4) 



Meyerhof showed in Ref. 7 that 



(cr 



J' 

ri' 



sin 25^, ^ 



( B ) « 0 



( 20 ) 



Also 




sin 26^, 
cos 26^1 



) s ( A ) «0 



on the same basis as the proof of (20). Therefore, A® = B® « 0 and 
expression ( 0 ^+) is reduced to 

< AAAa + ABAb ) (C5) 

where we have also used ( a®A A ) = ( b°A B ) = 0. 

J * 

A A and A B are functions of 5 i while A a and A b refer only 

Ju 

to 6 ^. For Z f Z' and/or J / J', the fluctuation over the averag- 
ing interval in functions of are upcorrelated to those in func- 

tions of thus, expression (C5) is equal to zero, and we have 

proved Eq. (Cl) for this case. 



For J = J' 



and 



the LHS of Eq. (Cl) reduces to 



< 




{ 



Z = Z' 

2 sin^ 
sin 2 6 




) , 



and we may use the 


result proved in Ref, 


7, 




^ . 2 J' 






JZ 


2 sin 6 , . 
1 • 






‘^ri’ 


^ sin 2 6^, ■ 


0 J 



which completes the proof of Eq. (Cl). 



APPENDIX D 



GAMMA RAY SUPPRESSION PREAMPLIFIER 

A space charge saturation gamma-ray suppression preamplifier 
vas used to distinguish neutron induced from gamma-ray induced pulses 
in the detector. The circuit is shovn in Fig. 22. It is basically 

32 

that of Owens, modified to include more external controls, a diode 
clipper on the nonlinear output, and an integrating network for the 
remaining positive overshoot. The external control on resistor R^^ 

33 

was siaggested by the work of Funsten and Cobb. 

Certain organic scintillators emit light under particle bom- 
bardment, which is composed of an initial spike of fluorescence with 
about sec decay time, followed by one or more slower components 

34 

of greater than 200 m\i sec decay times. Ion recombination may ac- 

35 

count for the slow components. The ratio of the amplitude of the 
initial spike to that of the first slow component is a function of the 
type of particle causing the fluorescence. 

The circuit in Fig. 22 is designed to produce linear voltage 
pulses, i.e., voltage pulses which are proportional to the amplitude 
of the initial fluorescence, at dynode. 10. The remainder of the re- 
sistor chain is arranged deliberately to produce space-charge satura- 
tion in the region between dynode Ik and the anode, during the initial 
spike. When the slow- component cascade electrons arrive in the region 
of space-charge saturation they are sufficiently repelled to produce a 
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Fig. 22. Circuit diagram of a space charge saturation gamma 
ray suppression preamplifier. The four externally controlled elements 
are marked "FOC," "ACC," "R^/' and "dYN ik." "UN" means linear 
pulse output. 'V supp" means the gamma ray suppression (nonlinear) 



output . 



126 



SPACE CHARGE SATURATION y SUPPRESSION PREAMPLIFIER 
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positive overshoot pulse at dynode ih. Thus^ the height of the positive 
overshoot is a measure of the ratio of slow component to initial spike 
amplitude and it is possible to use this effect to distinguish the type 
of particle which initiated the fluorescence. 

One is particularly interested in distinguishing neutron- 
initiated pulses (i.e.^ from recoil protons) from pulses caused by 
gamma-rays (i.e., due to electrons), and the method has proved useful 
for this purpose. 0\ir particular circuit is arranged such that the 
negative pulse at dynode l4, caused by the initial spike for either in- 
cident neutrons or gamma-rays, is clipped by a diode. The remaining 
positive overshoot is, of coiirse, no longer a linear effect, and we 
have chosen to integrate this pulse in an attempt to maximize the dis- 
tinction between neutron and gamma-ray initiated pulses. 

The circuit was effective in gamma suppression for neutron 
energies greater than about 0.3 Mev; however, some care must be taken 
in bias sing out gamma-ray induced pulses and to select neutron caused 
pulses, particularly when working at low neutron energies. 

An example of the nonlinear spectrum obtained at 1.7-Mev neu- 
tron energy is given in Fig. 23. A nonlinear spectrum of gamma- ray 
initiated pulses is also shown for comparison. These gamma- rays were 
produced under Van de Graaff machine operating conditions with the 
proton beam incident upon a stopper in front of the tritiimi target. 

While working at low neutron energies (<0.5 Mev) we found that 
very small dynode Ik voltages were necessary to provide good gamma-ray 
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Fig. 23. Nonlinear pulse spectra obtained at = I.7 Mev 

from a space charge saturation gamma ray suppression preamplifier. 

The detector crystal was in the 0 “ neutron beam direction. The 

solid circles show the nonlinear spectrum with the proton beam on 

the target (T(p,n)He'^) . Open circles show the nonlinear spectrum 

with the proton beam blocked by a tantalum beam stopper in front of 

3 

the T target. 
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suppression. This produced smaller nonlinear output pulses, which 
could he compensated by increasing the high voltage on the photomulti- 
plier, but which eventually represented an essential limitation on the 
circiiit fianction because of the increasing influence of tube noise. 



APPENDIX E 



MULTIPLE SCATTERING AND FINITE SCATTERER CORRECTIONS 

The slab scatterers used in the experiment were between 1/5 
and 1/2 neutron mean free paths in thickness. Thus a significant 
nmber of multiply scattered neutrons were detected. Cross section 
data must be corrected for this effect. A second correction was re- 
quired because the detection crystals were close enough to the slab 
scatterer to detect once-scattered neutrons over an angular range of 
about ±8 degrees in polar angle with respect to the line of crystal- - 
scatterer centers. One can see this by examining Fig. 1 which shows 
the scatterer and detectors in typical positions. Because the detectors 
were sampling the differential cross section over this angular spread, 
peaks in the experimental angular distribution were reduced and valleys 
filled. We will refer to the correction for this effect as a "finite 
scatterer correction." 

We were not prepared to make the computer analysis necessary 
for obtaining accurate corrections because for cusp analysis a first 
order estimate was considered satisfactory. 

The measurements of differential cross sections versus energy 
were made for only one or two angles at a time . Hence, the experimental 
points were correlated in energy but not in angle . The absolute error 
in the (uncorrected) measured cross section was approximately ±15 per- 
cent (see Error Analysis) and was independently established for each 
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angle. Therefore, when one cross plots to present the angular distri- 
bution for different energies, the points determining the shape of the 
angular distribution are separately uncertain by ±15 percent. 

The shape of the differential scattering cross section is, of 
coxirse, essential to a determination of the multiple scattering correc- 
tion. We used first order m\altiple scattering and finite scatterer 

corrections in this experiment in two ways. We corrected the optical 

17 

model prediction of Moldauer for the angular distribution of cerim 
and compared the results with experimental data (see Fig. 8). We also 
corrected the calculated "energy-averaged" cusp effects before applying 
them to the experimental data for Ce, Ba, and Zr. 

In the Fe experiments, the structure in the data did not justify 
any detailed estimate of the corrections, short of a full computer study 
on the angular and energy dependence of multiple scattering. On the 
basis of characteristic scatterer dimensions we estimated an upper 
limit of 20 percent for the corrections in the case of iron. 

In the lithium experiments, the scatterer was 0.17 neutron 
mean free paths in characteristic dimension. This dimension is small 
in comparison with other scatterer thicknesses typically used in this 
experiment. The angular distribution is peaked backwards and the lab 
neutron energy loss for back scattering is from 30 to 40 percent of 
incident energy (below the detection bias of the experiment). All of 
these factors, particularly the last, contribute to an estimate that 
the correction for multiple scattering in this case is small. The 
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general agreement in magnitude between the measured lithium angular 
distribution and other work is shown in Fig. 9 and supports this con- 
clusion . 

It remains to discuss in detail the first order multiple scat- 
tering and finite scatterer corrections for Ce, Ba, and Zr. We used 

36 

a procedure and notation presented by Nauta. Three assumptions 

is independent 

of K, the order of scattering, where is the probability of scat- 

tering elastically K times before leaving the scatterer; second^ the 
angular distribution of multiply scattered neutrons is the same as in 
an infinite medium of the same material; and thirds there is no loss 
of lab energy upon scattering. The scatterer is characterized by a 
parameter 6 which is the fraction of neutrons which^ having been once 
scattered elastically, have undergone at least one more interaction. 

Neglecting third and higher order scattering, the necessary 

^ 36 

formulas are : 

da = f (0) d 0 3 [1 + (S ®^(0)/s da (0) (El) 

corr m.s. exp 2 1 exp 

where 

00 

S ®\0)/s ®^(0) = (2K6/0 d 0 ) Z [2/(2i + l)](a.)^ P.(cos 0) 

J. ex X X 

1 1 
= 2(21 + 1 ) / d 0 ^^ 

and where d0 is the true differential cross section; d0 is 

corr ' exp 

the measured differential cross section; 0,j, is the total cross section; 



underlie the calculation; first, the ratio 
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and with =/ S^®^(0)dn. The last 

definition emphasizes the fact that these equations are circularly 

related and that approximations of some kind are necessary to obtain 

a solution. Our procedure was to substitute an optical model angular 

17 

distribution for thus we calculated a correction to the opti- 

cal model prediction for comparison with experiment. Satisfactory . 
agreement between the corrected optical model distributions and experi- 
ment was considered justification for applying these corrections to 
the comparison of predicted cusp effects with experiment. This general 
procedure is analogous to that usually adopted in solving Eq. (El); 
one assimes a function ^7 iteration of Eq. (El) forces it 

to converge to ^^^orr* stopped this process at a first iteration. 

In essence, the method separates that part of the multiple scat- 
tering which depends upon the angular distribution characteristic of 
the scattering material from a part which depends upon the size, shape, 
and aspect (with respect to the beam axis) of the scatterer. The 
latter dependence is contained in the parameter 5. 

Lacking a better estimate of 5, we assumed that the quantity 

9l/<Pii in the detailed work of A. Langsdorf, Jr., R. 0. Lane, and 

37 

J. E. Monahan is 1-8, where is the emergent singly scattered 

fl\ix per unit solid angle and is the unit idealized scattered 

flux per unit solid angle. The parameter 5 is dependent upon the 
angle of scattering at which the cross section is measured because 
slab scatterers are extremely asymmetric . In the work of Nauta , ^ who 



used ring scatterers, 8 could be assumed independent cf angle. 
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Accordingly, we used the values of given in Fig. 9 of 

Ref. 37 lo obtain a 6 for each scatterer and for each angle of inter- 
est. In most cases extrapolation of the curves was necessary. These 
cxirves are entered with a particular value of h, the scatterer thick- 
ness, in units of neutron mean free path. 

Our scatterer thicknesses, calculated at the neutron energies 
of interest, were'; for Zr, H = 0.32; for Ce, h = 0.U6; and for Ba, 
h = 0 . 55 . 

The moments a^ of the optical model angular distribution were 
included to a|^. Multiple scattering corrections as a function of 
angle were then computed using Eq. (El) and the appropriate values of 
B . The degree to which corrected optical model predictions agreed 
with experimental data is indicated in Fig. 8. We estimated the over- 
all uncertainty in the absolute multiple scattering corrections result- 
ing from the assumptions and approximations in the foregoing treatment 
to be ±30 percent. Expressed as the fraction g ^ which is con- 
venient for modifying data, the vincertainty in multiple scattering 
corrections adds to the uncertainty obtained in the experimental error 
analysis. The uncertainty in f^ ^ varies with scatterer and detec- 
tion angle, as indicated in Table III. 

We nov consider the calculation of finite scatterer corrections . 
At each mean polar angle, 0^, once scattered neutrons are detected 
over an angle spread of 0^ ± A 0. The magnitude of 2 A 0 was de- 
temined graphically from the geometry of the experiment, see Fig. 1, 
and was found to vary from 10 to 35 degrees. The angular spread is a 
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function of 0 because the scatterer-to-detector distance varied 
o 

depending upon the experimental situation (counting rate, background to 

effect ratio) and also because the projected area of the scatterer, 

normal to the scatterer-detector line of centers, varies with 0 . 

o 

As in the correction for multiple scattering, the procedure 
was to correct an optical model prediction for experimental effects. 

The finite scatterer correction is defined by 



da 

corr 



= f„ (0) da 
f.s. exp 



(E2) 



f (0) was computed using 

X • S • 

N 

(®) * da /(1/N Z da .(sin 0 /sin 0 ) (E3) 

f.s. ' opt' ' opt' n' o 



where optical model prediction at mean angle 0^, 

and the denominator represents an average of the optical model distri- 
bution over N equal divisions of the spread 2 A 0 weighted by the 
solid angle subtended by the nth angle division. N was selected to 
give a representative average and varied from 5 to 8 depending upon 
the size of 2 A 0. 

The resxilts of this calculation are given in Table III. The 
magnitude of typical finite scatterer corrections is indicated on the 
angular distribution shown in Fig. 8, where it can be seen that the 
correction is small except in the valley near 0^ = 90 '^* 
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TABLE III. Multiple Scattering and Finite Scatterer Corrections 





e 

o 


f 

m.s . 


.s. 


Uncertainty 

in f X 

m.s. f .s . 

(percent) 


Cerium 


30“ 


0.95 


l.o4 


3 




50° 


0.87 


0.98 


5 




70° 


0.80 


0.92 


8 




9.0° 


0.75 


0.90 


10 




110“ 


0.78 


0.99 


9 




130“ 


0.80 


1.00 


7 




150° 


0.85 


1.05 


6 


Barium 


90“ 


0.65 


0.90 


18 




125° 


0.70 


1.00 


l4 


Zirconium 


60° 


0.90 


~1.0 


k 




130“ 


0.85 


~1.0 


6 


Iron 


All 


0. 8-1.0 


~1.0 


~10 


Lithium 


All 


~1.0 


1.0 


~ 5 


The 


factors in this 


table may 


be applied to 


the' .;experimental 



cross section data presented in this report, by using 



da = f f 
corr m.s. 
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